AO-A047  504 


UNCLAbSlFlEO 


GENERAL  ELECTRIC  CORPORATE  RESEARCH  AND  DEVELOPMENT  —ETC  F/G  11/2 
INVESTlbATlON  OF  CERAMICS  FOR  HIGH  TEMPERATURE  TURBINE  COMPONEN — ETC(U) 
JUN  77  C A JOHNSON.  S PR0CHA2KA  N62269-76-C-0243 


SRD-77-145 


NADC-75228-30 


ADA047504 


NADC- 75228-30 


H 

Qu- 

O 

O 


§ 


INVESTIGATION  OF  CERAMICS 
FOR  HIGH  TEMPERATURE  TURBINE  COMPONENTS 


C.  A.  Johnson  and  S.  Prochazka 
GENERAL  ELECTRIC  COMPANY 
Corporate  Research  and  Development 
P.  O.  Box  8 

Schenectady,  New  York  12301 


June  1977 


Final  Report  for  Period  25  March  1976  - 25  March  1977 


APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


Prepared  under  Contract  N62269-76-C-0243 


NAVAL  AIR  DEVELOPMENT  CENTER 
Warminster,  PA  18974 


for 

NAVAL  AIR  SYSTEMS  COMMAND 
Department  of  the  Navy 
Washington,  DC  20361 


Vt)  D C 

rWE®EISl[IlZfR\ 
UJ  DEC  U tsn  I 11 


JSEDTnilii: 


SRD-77-14S 


UNCLASSIFIED 

^>eteA»SIFICATION  OF  THIS  PACE  (Whan  P«t«  Bnfrmt/) 

I c/  ) REPORT  DOCUMENTATION  PAGE 


l2.  GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 
S.  RECIPIENT'S  CATALOG  NUMBER 


rNADCf|75228-3j 


4.  TITLE.r«n4tJluUillR) — 

JNVESTIGATION  OF.  CERAMICS  FOR  inGH 
TEMPERATURE  TUiTbINE  COMPONENTS  * 


jO)  Ciirtis  A./johnson  / 

Svante/Prochazka  / ^ 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

general  ElectriCj^OompRHy 
v Corporate  Research  and  Development 

S:;henectady,  NY  123oi ^ 

n.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  7TT 

Dept,  of  the  Navy  ( /' 

Naval  Air  Systems  Command 
Washington,  DC  20361 

14.  MONITORING  AGENCY  NAME  • AODRESSflf  (»/f«r«nl  fraai  CcntrolUni  OUlet) 

Dept,  of  the  Navy 

Naval  Air  Development  Center 

Warminster,  PA  18974 


S.  TYPE  OP  REPORT  9 PERIOD  COVERED 

Final  Report  - 

76  Mar  25  to  77  Mar  25 

k.  pgBFnniiiiir  nntt,..wpf»R-i‘  uiiutp 

/ SRD-77-145  ( 

“Br^^lWTtPTWrORANT  NUM8ERC*; 

^N62269-76-C-.<)243 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 

A03V3200/001B/6F54592201 


13.  NUMBER  OF  PAGES^~-“  - 

Tr""iECURTfY'CLASirfo77Si7«porir 

Unclassified 


15«.  OECLASSIFICATION/DOWNGRADING 
SCHEDULE 


I IS.  distribution  statement  fo/  thl9  Report) 


Approved  for  Public  Release;  Distributi<m\Unlimited 


17.  distribution  stj 


18.  SUPPLEMENTARY  NOTES 


1 19.  key  words  (Contintf  on  1/  n«c««««ry  and  id^nUfy  by  block  number) 


silicon  carbide,  sintering,  fracture  strength,  crack  blunting,  grain  growth, 
oxidation,  fracture  statistics,  Weibull  statistics,  fracture  origins,  stress 
rupture,  statis  fatigue,  dynamic  fatigue,  creep,  fracture  toughness 


2o!  a^^RACT  (Contlnu0  on  »ldm  It  n0C9»aary  and  Identity  by  Mock  nttmbor) 

The  mechanical  properties  of  boron-doped  sintered  SiC  were  studied  as  a 
function  of  several  fabrication  and  testing  variables.  The  reduction  in  severity 
of  surface  defects  due  to  the  dissociation  of  SiC  at  high  temperatures  was  charac- 
terized as  a function  of  furnace  pressure  during  sintering  and  as  a function  of 
refiring  temperature  after  machining.  The  to polytype  transformation  and 
the  accompanying  grain  growth  were  found  t^decrease  the  strength  & maximum^ 

(Continued  on  bpck)  At-.p^A 


FORM 
I JAN  71 


EDITION  OF  1 NOV  88  IS  OBSOLETE 


UNCLASSIFIED 


ml  TMi  ilfttB  pbib  IribuB) 


SCquWITY  CLASSiriCATlOH  OF  THIS  PAOEOrhan  Data  Entarad; 


of  approximately  25  percent.  Oxidation  for  times  to  2000  hours  at  150^^  in 
air  resulted  in  no  net  loss  of  strength.  The  time -tempera^re -strength  inter 
dependence  was  studied  at  temperatures  from  1200  to  2100^  and  times  to 
lailure  from  one  seco^  to  five  hours.  The  creep  behavior  was  also  studied 
in  the  vicinity  of  210(5^.  A method  was  demonstrated  to  derive  statistical 
quantities  such  as  the  Weibull  modulus  from  the  positions  of  fracture  origins, 
and  two  examples  of  fracture  data  were  discussed  which  illustrate  the  effects 
of  bimodal  flaw  distributions.  Hot  pressing  experiments  suggest  that  the 
large  cavities  often  observed  in  sintered  microstructures  are  caused  by 
bridging  of  ^^lomerates  during  powder  compaction.  Preliminary  studies  of 
increasing through  the  addition  of  second  phase  particles  has  not  been 
successful.^^ 


fC  S'va 


UNCLASSIFIED 


TABLE  OF  CONTENTS 


I Introduction  1 

II  Effect  of  Sintering  Furnace  Pressure  on  Fracture  Strength.  ...  6 

III  Effect  of  Refiring  after  Machining  on  Fracture  Strength 9 

IV  Effect  of  3 to  ot  Transformation  on  Fracture  Strength 12 

V Effect  of  Oxidation  on  Fracture  Strength  20 

VI  Fracture  Statistics 29 

A.  Derivation  of  Statistical  Parameters  of 
Fracture  from  the  Locations  of  Fracture 

Origins 35 

B.  Multiple  Flaw  Distribution 48 

VII  Time-Temperature  Dependence  of  Fracture  Strength  61 

VIII  Creep  Measurements  73 

IX  Origin  of  Strength  Limiting  Flaws 88 

X Increasing  Fracture  Toughness  Through  the  Addition 

of  Second  Phases lOO 

References 

Appendix  A 

Appendix  B 

t 


iii 


FOREWORD 


This  work  was  sponsored  by  the  Naval  Air  Development  Center  and 
carried  out  in  the  General  Electric  Corporate  Research  and  Development 
Center,  Physical  Chemistry  Laboratory,  Ceramics  Branch  - Dr.  R.  J. 
Charles,  Manager,  under  Contract  No.  N62269-76-C-0243.  It  is  a 
continuation  of  an  investigation  in  high  temperature,  high  strength 
ceramics  for  structural  engineering  application  which  was  initiated 
in  1971.  This  portion  was  performed  during  the  period  from  March  25, 
1976  to  March  25,  1977.  Mr.  I.  Machlin  of  the  Naval  Air  Systems 
Command  was  the  program  consultant. 

The  authors  would  like  to  acknowledge  the  valuable  technical 
discussions  with  Drs.  R.  A.  Giddings,  R.  J.  Charles  and  G.  G. 

Trantina.  The  technical  assistance  of  C.  Bobik  is  also  appreciated. 
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SUMMARY  OF  IMPORTANT  RESULTS 

1)  Bend  specimens  of  SiC  sintered  at  low  furnace  pressures  (0.3  to  0.4 
torr)  were  found  to  have  average  strengths  which  were  25  to  30  percent 
higher  than  similar  specimens  sintered  under  one  atmosphere  of  inert 
gas.  The  increase  in  strength  at  low  furnace  pressures  is  believed  to 
result  from  a crack  "blunting"  mechanism  through  the  dissociation  of 
Sic  at  high  temperatures  and  low  pressures. 

2)  Sintered  and  machined  specimens  of  SiC  were  found  to  increase  in 
strength  as  much  as  20  percent  by  refiring  them  near  the  original 
sintering  temperature  at  a furnace  pressure  of  0.3  to  0.4  torr.  This 
effect  is  also  believed  to  be  the  result  of  a crack  blunting  mechanism. 

3)  The  bend  strength  of  sintered  SiC  was  found  to  be  a relatively  insensitive 

( 

function  of  the  degree  of  3 to  a transformation  despite  an  accon^janying 
increase  in  grain  size.  The  average  strength  of  fully  transformed  a-SiC 
was  approximately  25  percent  lower  than  that  of  a fine-grained  3-SiC 
I microstructure . 

4)  Specimens  of  sintered  SiC  which  were  oxidized  in  still  air  at  1500°C 
for  times  up  to  2000  hours  showed  a nearly  parabolic  rate  of  oxidation 
and  no  net  loss  of  strength. 

5)  A method  of  quantitatively  determining  the  statistical  parameters  of 
fracture  from  the  positions  of  fracture  origins  was  derived  for  two- 
parameter  Weibull  statistics  and  demonstrated  on  a series  of  abraded 
glass  rods. 

6)  Two  sets  of  literature  data  were  reanalyzed  to  quantitatively  illustrate 

examples  of  multiple  flaw  distributions.  In  both  cases,  the  observations 
were  more  accurately  described  by  assuming  two  active  flaw  distributions 

than  the  more  conventional  assun^tion  of  only  one  flaw  distribution. 
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7)  The  stress  rupture  behavior  of  sintered  SiC  was  studied  in  biaxial 
tension  at  temperatures  of  approximately  2100®C  by  pressurizing 
slip  cast  tubes  with  nitrogen.  The  results  were  correlated  with 
previous  measurements  of  stress  rupture  in  bending  at  lower 
temperatures . 

8)  The  time-temperature  dependence  of  fracture  strength  in  sintered 
SiC  was  also  studied  in  bending  as  a function  of  temperature  to 
IbOO^C  and  as  a function  of  loading  rate  at  1600®C.  At  a total 
time  to  failure  of  five  minutes,  the  strength  was  found  to  be 
independent  of  temperature  to  approximately  1200'^C  and  dropped  to 
approximately  85  percent  of  the  room  ten^jerature  strength  by  1600®C. 
Quantitative  analysis  of  strength  degradation  due  to  subcritical 

crack  growth  at  1600®C  was  conqjlicated  by  concurrent  surface  degradation 
due  to  oxidation. 


9)  Hie  creep  rate  of  sintered  SiC  was  measured  at  temperatures  in  the 
vicinity  of  2100'’C  using  the  pressurized  tubes  mentioned  above  and 
resulted  in  an  activation  energy  for  creep  of  246  Kcal/mole.  By 
assuming  that  the  primary  mechanism  of  material  transport  was  volume 
diffusion,  these  data  along  with  previous  measurements  of  creep  in 
bending  at  1600®C  were  used  to  derive  an  apparent  diffusion  coefficient 

10 

D^pp  = 3.25  X lO^'^  exp  [-(225  + 29) /RT] 

10)  A careful  comparison  of  the  properties  of  boron-doped  hot-pressed  and 
sintered  SiC's  has  led  to  the  conclusion  that  the  large  voids  often 
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observed  in  sintered  microstructures  are  caused  by  bridging  of 
agglomerates  during  the  compaction  step  of  specimen  fabrication. 

11)  Preliminary  studies  of  increasing  the  fracture  toughness  of  sintered 
SiC  through  the  addition  of  particulate  second  phase  dispersions  have 
identified  several  second  phase  systems  which  are  compatible  with  the 
sintering  process.  The  only  measurements  of  on  two  phased  materials 
to  date,  however,  revealed  no  distinguishable  change  in  toughness  with 
the  addition  of  S w/o  of  B-C. 
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Major  steps  in  the  improvement  of  the  gas  turbine  engine  have 
been  associated  with  advancements  in  superalloys  and  cooling  technology. 
Although  further  improvements  are  expected,  definite  performance 
limitations  exist  in  traditional  metallic  systems,  and  a break- 
through requires  substantially  increased  component  operating  temperatures. 
The  potential  of  ceramics  for  such  a breakthrough  has  been  recognized 
and  a considerable  effort  has  been  launched  to  develop  and  evaluate 
new  high  temperature  structural  ceramics  as  well  as  to  develop  the 
design  technology  necessary  for  their  engineering  application. 

Silicon  carbide  and  silicon  nitride  have  been  identified  as 
promis'  andidates  based  on  properties  such  as  oxidation  stability, 
ere  ance,  thermal  shock  properties,  strength,  availability 

a r criteria.  A number  of  consolidation  techniques  have  been 

developed  for  SiC  and  Si^N^  and  have  resulted  in  a variety  of 
materials.  Such  diversification  is  highly  desirable  since  one 
material  cannot  satisfy  all  of  the  emerging  requirements. 

The  objective  of  this  continuing  program  has  been  to  develop 
the  performance  characteristics  of  a fabricable  form  of  silicon 
carbide  such  that  an  operational  temperature  approaching  3000®F  may 
be  realized.  The  critical  aspects  of  performance  include  strength, 
variability  of  strength,  stress  rupture  behavior,  fracture  toughness 
and  chemical  stability  in  oxidizing  environments.  It  must  be 
anticipated  that  the  precise  performance  levels  required  of  the 
material  in  these  categories  are  a strong  function  of  the  design 
innovations  that  can  be  achieved  for  nonductile  materials.  In 
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general,  however,  the  materials  development  must  maintain  the  ’ ’ 

objectives  of  developing  and  maintaining  very  high  temperature 

i 

strength  and  surface  integrity  without  sacrificing  fabrication 

I 

capability  or  creep  resistance. 

In  the  first  year  of  this  program,  the  system  SiC-B  was 

identified  as  a promising  candidate  for  the  development  of  good 

oxidation  resistance  and  strength  at  high  temperatures . The  work 

was  carried  out  under  Contract  No.  N00019-71-C-0290  and  discussed 

in  detail  in  report  SRD-72-035. 

During  the  second  year  it  was  shown  that  the  elimination  of 

either  silicon  or  silica  from  hot-pressed  SiC  through  the  addition 

of  carbon  resulted  in  significant  improvements  of  high  temperature 

strength,  creep  resistance,  and  stress  rupture  properties.  Bend 

strengths  of  about  75  kpsi  at  1550®C  and  creep  rates  at  the  same 

-9 

temperature  of  less  than  10  /sec  at  25  kpsi  were  achieved.  The 
potential  of  a boron-doped  single  phase  SiC  ceramic  as  a high 
temperature  structural  material  was  firmly  established.  This 

i 

second  period  of  the  program  was  supported  by  Contract  No.  N00019-72-  * 

C-0129  and  reported  in  SRD-72-171.  | 

Work  in  the  third  year  of  the  program  was  carried  out  under  | 

Contract  No.  N62269-73-C-0356  and  was  reported  in  SRD-74-040  of  | 

April  1974.  It  was  demonstrated  that  conventional  sintering  of  | 

Sic  powders  to  high  density  was  possible  and  that  mechanical  pro- 
perties identical  to  those  of  the  earlier  hot-pressed  materials  could 
be  obtained.  In  addition,  it  was  shown  that  the  new  sintering  process 

i 

was  amenable  to  fabrication  of  complex  parts  by  several  conventional  ^ j 

fabrication  techniques . 
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The  work  during  the  following  year  was  carried  out  under 
Contract  No.  N62269-74-C-0255  and  summarized  in  report  No.  SRD-74- 
123.  The  objectives  included:  a study  of  the  effect  of  additives, 
impurities  and  sintering  conditions;  the  development  of  fabrication 
capability  of  SiC  on  a laboratory  scale;  and  the  evaluation  of  the 
physical  and  mechanical  properties  of  the  improved  SiC.  A more 
detail'*d  investigation  of  the  P to  a transformation  in  SiC  was 
necessary  to  obtain  control  over  exaggerated  grain  growth  during 
sintering.  Thermal  expansion,  elastic  modulus,  thermal  conductivity, 
modulus  of  rupture,  Weibull  modulus,  fracture  energy  and  creep  were 
evaluated  on  specimens  prepared  in  the  laboratory. 

During  the  fifth  year,  the  program  was  supported  under  Contract 
No.  N62269-75-C-0122  and  reported  in  SRD-75-126.  Work  carried  out  in- 
cluded the  development  of  improved  powder  processing  techniques  where 
the  SiC  powder  was  dispersed  in  a liquid  such  as  benzene  and  passed 
through  a very  fine  screen  to  remove  large  agglomerates  and  foreign 
material.  The  liquid  was  then  removed  by  spray  freezing  and  freeze 
drying  to  minimize  the  segregation  of  soluble  binders  and  lubricants. 

The  dependence  of  density  and  microstructure  on  sintering  variables 
such  as  temperature,  time,  green  density,  particle  size  and  sintering 
atmosphere  were  studied  quantitatively  in  a high  temperature  dialatometer . 
The  fracture  strength  was  studied  as  a function  of  fabrication  and 
testing  variables  such  as  powder  processing  technique,  density, 
sintering  atmosphere,  testing  temperature,  etc.  Testing  was  also 
carried  out  under  conditions  of  thermally  induced  stresses  to  study 
thermal  shock  and  thermal  fatigue  of  sintered  SiC. 
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During  the  present  contract,  efforts  have  been  directed  toward 
three  general  areas:  the  continued  characterization  of  those  properties 
necessary  for  the  effective  design  of  high  temperature  structural  components 
from  sintered  SiC;  the  evaluation  and  control  of  microstructural  inhomo- 
geneoties  which  are  commonly  associated  with  fracture;  and  the  initial 
feasibility  study  of  increasing  the  fracture  toughness  of  SiC  through 
the  addition  of  appropriate  second  phase  particles.  Since  sintered  SiC 
is  being  developed  as  a structural  material,  it  is  appropriate  that  all 
three  of  these  areas  of  study  are  concerned  with  the  characterization, 
understanding  and  optimization  of  strength  related  properties.  Of  the 
three  areas,  the  first  has  been  emphasized  during  the  previous  year.  The 
fracture  strength  was  studied  as  a function  of  several  fabrication  and 
testing  variables  which  must  be  optimized  to  assure  the  highest  pro- 
bability of  success  and  must  be  sufficiently  well  characterized  to 
satisfactorily  predict  the  behavior  of  a given  component  design  under 
service  conditions . The  improvement  in  strength  caused  by  blunting  of 
surface  flaws  through  the  dissociation  of  SiC  at  high  temperatures  was 
characterized  both  during  the  initial  sintering  operation  where  the 
severity  of  surface  fabrication  defects  is  reduced  and  during  a refiring 
operation  after  machining  where  the  severity  of  machining  induced  surface 
damage  is  also  reduced.  The  fracture  strength  was  characterized  as  a 
function  of  the  degree  of  6 to  a transformation  and  as  a function  of 
exposure  time  in  air  at  1500°C.  The  statistics  of  fracture  were 
approached  by  developing  and  verifying  the  tools  to  satisfactorily 
describe  the  size  dependence  of  fracture  strength  in  sintered  SiC  as 
well  as  in  other  potential  structural  ceramics.  The  temperature-time- 
strength  interdependence  was  studied  in  sintered  SiC  from  1200  to  2100°C 
and  at  times  to  failure  of  one  second  to  five  hours.  In  the  area  of 
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the  evaluation  and  control  of  microstructural  inhomogeneoties , 
experiments  have  supported  the  hypothesis  that  the  large  voids 
which  commonly  serve  as  strength  limiting  flaws  are  caused  by 
massive  bridging  of  weak  agglomerates  during  the  powder  pressing 
step  of  specimen  fabrication.  Finally,  in  the  area  of  fracture 
toughness,  second  phase  particles  have  been  incorporated  into 
sinterable  SiC  powder,  the  mixtures  have  been  hot-pressed  or 
sintered  to  high  density  and  the  resulting  materials  have  been  tested 
at  room  temperature. 

A concurrent  program  is  oriented  toward  the  synthesis  and 
characterization  of  sinterable  silicon  carbide  powders.  These  results 
are  not  reported  here  but  will  become  available  elsewhere.  In  order 
to  maintain  continuity  with  previous  reports,  the  majority  of  work 
reported  below  was  carried  out  with  powders  from  Pittsburgh  Plate 
Glass  Company  which  were  prepared  by  a proprietary  process  from 
chlorinated  silicon  compounds  and  hydrocarbons.  The  characteristics 
of  these  powders  have  been  reported  in  Ref.  1,  p.  62. 
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II.  EFFECT  OF  SINTERING  FURNACE  PRESSURE  ON  FRACTURE  STRENGTH 


r 
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Included  in  last  year's  final  report  were  the  results  of  a 
limited  study  of  the  effect  of  furnace  atmosphere  pressure  during 
sintering  on  the  fracture  strength  of  bend  specimens  tested  in  the 
"as-sintered"  state.  These  data  indicated  that  higher  strengths 
could  be  obtained  by  sintering  in  vacuum  (.2-. 4 torr)  than  by  sin- 
tering in  one  atmosphere  of  inert  gas,  despite  equivalent  grain 
sizes  and  densities.  Other  experiments  described  in  the  same 
reference  indicated  that  the  origin  of  the  higher  strength  of  speci- 
mens sintered  at  low  pressure  was  the  "blunting"  of  cracks  and 
pressing  flaws  through  the  dissociation  of  SiC  from  all  free  surfaces. 

The  resulting  silicon  vapor  is  believed  to  travel  to  cooler  parts  of 
the  furnace  where  it  recombines  with  carbon  insulation,  etc.  The 
carbon  remains  on  the  original  surface  and  is  readily  observable  on 
all  specimens  sintered  at  low  pressure.  The  kinetics  of  the  dissoci- 
ation are  strongly  influenced  by  the  total  pressure  of  inert  gaseous 
species  above  the  surface  of  the  reaction.  At  low  furnace  pressure  the 
reaction  proceeds  more  rapidly  than  at  higher  pressure,  thus  dissociating 
a greater  quantity  of  SiC  and  more  effectively  "blunting"  any  surface 
cracks,  etc. 

This  phenomenon  has  been  judged  to  be  sufficiently  important  to 
warrant  a more  careful  characterization  of  the  strength  as  a function 
of  furnace  pressure  during  sintering.  A series  of  150  specimens  were 
die-pressed  in  a 0.2  x 2.0  in.  die  at  four  kpsi  (28  MPa)  using  boron- 
doped,  B-SiC  powder  which  had  been  screened  through  an  8 pm  screen  as 
a dispersion  in  benzene,  spray  frozen,  and  freeze  dried  according  to 
the  procedures  described  previously  . The  bars  were  sintered  in  ten 


groups  of  fifteen  specimens  in  a carbon  resistance  furnace  at 
pressures  ranging  from  0.35  to  760  torr  of  helium.  Pressure  was 
controlled  with  a cartesian  diver  pressure  control  with  a small, 
continuous  flow  of  helium.  All  specimens  were  heated  to  temperature 
at  1000®C  per  hour  and  held  at  2085 ®C  for  30  minutes,  resulting  in 
densities  from  94  to  95  percent  based  on  a theoretical  density  of 

3 

3.21  g/cm  . They  were  tested  in  the  "as-sintered"  state  (no  machining, 

etc.)  in  three-point  bending  on  a one  and  one-half  inch  span  at 

room  temperature  using  the  sphere-loaded  test  fixture  described 
2 

previously.  The  results  are  plotted  on  Figure  1 as  fracture 
strength  vs.  furnace  pressure.  Each  data  point  represents  the 
average  strength  of  one  group  and  the  error  bars  include  the  95 
percent  confidence  limits  on  the  average  as  calculated  by  Student's 
t-distribution.  The  data  points  of  the  highest  and  lowest  furnace 
pressures  represent  two  groups  (30  specimens)  each. 

As  shown  in  Figure  1,  the  lowest  furnace  pressures  yielded  an 
average  strength  of  approximately  105  kpsi  (720  MPa)  and  the  highest 
pressure  yielded  approximately  80  kpsi  (550  MPa).  The  sigmoidal 
shape  of  the  dashed  curve  drawn  through  the  data  reflects  the  expected 
limiting  strengths  at  both  high  and  low  pressures.  As  furnace  pressure 
is  decreased,  the  "blunting"  of  surface  flaws  by  material  removal  can 
only  increase  the  fracture  strength  to  the  point  where  failure  is  con- 
trolled by  internal  flaws.  (Internal  defects  should  not  be  affected 
appreciably  by  changing  furnace  pressure.  Very  little  silicon  vapor 
can  escape  to  the  specimen  surface,  therefore  the  internal  partial 
pressure  of  silicon  is  always  near  equilibrium  with  effectively  zero 
rate  of  dissociation) . At  furnace  pressures  below  that  point  where 


failure  is  controlled  by  internal  flaws,  the  strength  is  expected 
to  be  independent  of  furnace  pressure  during  sintering.  The  location 
of  fracture  origins  supports  this  statement  with  an  increasing  pro- 
portion of  failures  observed  from  subsurface  flaws  as  the  furnace 
pressure  was  decreased.  A lower  limiting  strength  might  also  be 
expected  at  high  furnace  pressures  where  the  kinetics  of  the  dissoci- 
ation reaction  are  reduced  to  the  point  of  no  material  removal . ' This 
is  also  supported  qualitatively  by  the  observation  of  no  free  carbon 
on  the  surfaces  of  bodies  sintered  in  one  atmosphere  of  helium,  whereas 
similar  firing  conditions  at  low  pressure  produce  an  obvious  coating  of 
carbon  on  all  surfaces. 


PRESSURE  DURING  SINTERING  (TORR  OF  HELIUM) 

Figure  1.  Effect  of  Furnace  Pressure  During  Sintering 
on  the  Fracture  Strength  of  Specimens  Tested 
in  the  "As -Sintered"  State 
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III.  EFFECT  OF  REFIRING  AFTER  MACHINING  ON  FRACTURE  STRENGTH 

Another  phenomenon  first  described  in  last  year's  final  report 
was  the  use  of  vacuum  refiring  to  increase  the  fracture  strength  of 
machined  bodies.  Specimens;  which  are  dieimond  ground  are  often  observed 
to  fail  from  machining-induced  flaws.  Since  diamond  grinding  will  be 
unavoidable  for  many  structures  which  require  critical  dimensional 
tolerances,  a means  of  repairing  a large  degree  of  the  machining  damage 
without  disturbing  the  dimensions  will  be  very  useful. 

The  procedure  as  described  previously  involved  refiring  at  or  i 

I 

near  the  original  sintering  temperature  in  vacuum  (0.2-0. 4 torr)  for  j 

i 

20  to  30  minutes  and  resulted  in  a 30  to  40  percent  increase  in  the  j 

i 

fracture  strength  for  the  machined  specimens  reported  at  that  time.  i 

Refiring  at  or  very  near  the  original  sintering  temperature,  however,  j 

risks  both  a small  change  in  dimen.sions  due  to  additional  densification  | 

and  an  increase  in  the  largest  grain  size  due  to  exaggerated  grain  growth  I 

f 

during  any  g to  o transformation.  The  mechanism  of  strength  j 

i 

i 

improvement  by  this  refiring  technique  is  believed  to  be  the  same  dis-  i 

sociation  reaction  and  crack  "blunting"  mecheinism  discussed  in  the 

previous  section.  Since  the  dissociation  of  SiC  proceeds  more  slowly 

at  lower  temperatures,  however,  there  was  an  obvious  interest  in  | 

determining  the  lowest  refiring  temperature  which  could  be  used  to  ] 

repair  machining  damage.  It  was  therefore  decided  to  study  the  refiring  | 

of  machined  specimens  at  ten^jeratures  below  the  original  sintering 

temperature  in  the  hope  that  machining  damage  can  be  substantially 

repaired  without  the  risks  mentioned  above. 

A series  of  98  specimens  of  SiC  were  prepared  in  the  same  manner 
as  those  described  in  the  previous  section.  They  were  all  sintered  in 
vacuum  (.3-. 4 torr  of  helium)  at  approximately  2080“C  for  20  minutes 

d 
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and  were  diamond  ground  (220  grit)  from  approximately  0.165  in. 
sqiiare  cross  section  to  0.150  in.  square  cross  section.  They  were 
refired  in  vacuum  in  groups  of  12  at  six  different  temperatures  with 
the  remaining  26  tested  without  refiring  to  estimate  the  "as-machined" 
strength.  They  were  tested  as  described  in  the  previous  section  and 
the  results  are  plotted  on  Figure  2 as  fracture  strength  vs.  refiring 
temperature.  The  specimens  receiving  no  refiring  are  plotted  as  specimens 
refired  at  room  temperature.  Again  each  data  point  represents  the 
average  strength  and  the  error  bar  includes  the  95  percent  confidence 
limits  on  the  average.  The  strength  of  the  as-machined  specimens  averaged 
approximately  81  kpsi  (560  MPa)  and,  as  expected,  was  the  lowest  strength. 
The  highest  average  of  approximately  96  kpsi  (660  MPa)  was  measured  for 
those  specimens  refired  at  the  original  sintering  temperature. 

Despite  the  rather  large  degree  of  scatter,  a dashed  curve  was  drawn 
through  the  data  of  Figure  2 as  a suggested  relatxonsnip  .c*  ’een  strength 
and  refiring  temperature.  Again,  the  sigmoidal  shape  reflects  the  limiting 
effects  which  are  expected  at  very  high  and  very  low  refiring  temperatures. 
At  progressively  higher  temperatures  the  strength  will  become  limited  by 
internal  flaws  while  at  progressively  lower  temperatures  there  is  negligible 
dissociation  of  SiC. 

The  data  from  these  experiments  suggest  that  the  greatest  degree  of 
strength  recovery  in  machined  specimens  requires  a refiring  temperature 
no  lower  than  the  original  sintering  temperature.  It  is  expected,  however, 
that  the  increased  benefits  of  refiring  at  100°C  or  more  below  the 
sintering  temperature  will  more  than  compensate  for  the  small  sacrifice 
in  strength.  Another  variable  which  may  be  useful  but  has  not  been  charac- 
terized is  the  effect  of  longer  refiring  times  than  the  20  minutes  of  this 
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study.  Based  upon  the  present  understanding  of  the  phenomenon,  one  can 
safely  predict  that  longer  refiring  times  will  shift  the  entire  curve 
of  Figure  2 slightly  to  the  left  (comparable  strengths  achieved  at  lower 
temperatures  and  longer  times) . 


0.35  TORR  HELIUM 
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Figure  2,  Fracture  Strength  as  a Function  of  Refiring  Temperature 
for  Machined  Specimens 
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IV.  EFFECT  OF  3 TO  a TRANSFORMATION  ON  FRACTURE  STRENGTH 


r 

(i 


i 


The  cubic  polytype  of  SiC,  P or  3C,  is  generally  observed  to 

be  thermodynamically  unstable  relative  to  several  of  the  non-cubic 

a polytypes  in  the  temperature  regime  of  1800  to  2200‘’C.  The  phase 

transformation  which  tends  to  proceed  during  sintering  of  boron-doped 

3-SiC  powders  is  found  to  be  nucleation  controlled.  In  the  presence 

of  a small  number  of  favorable  sites  for  nucleation,  the  transformation 

can  lead  to  the  growth  of  extremely  large  crystals  of  a-SiC  which  often 

approach  1mm  in  size.  The  kinetics  and  microstructures  of  this  trans- 

1-3 

formation  have  been  reported  previously,  but  the  degradation  in 
strength  due  to  the  presence  of  such  large  grains  has  not  been  studied 
quantitatively . 

In  a study  of  fracture  strength  as  a function  of  degree  of  trans- 
formation during  sintering,  there  are  three  states  of  surface  preparation 
which  are  of  interest: 

• as-sintered 

• sintered  and  machined 

• sintered,  machined,  and  refired. 

Unfortunately,  it  is  difficult  to  prepare  specimens  in  the  as-sintered 

state  with  a uniform  degree  of  transformation  through  the  specimen 

thickness.  The  cross  section  of  a sintered  body  which  is  partially 

transformed  generally  consists  of  a thin  "skin"  .5mm)  of  primarily 

B-SiC  on  all  free  surfaces  with  the  remainder  of  the  body  uniformly 
2 

transformed.  The  most  viable  means  of  measuring  the  fracture  strength 
is  in  bending  where  the  highest  tensile  stresses  are  located  at  and 
near  the  surface  of  the  specimen.  The  results  of  bend  testing  in 
the  as-sintered  state  would  not  represent  the  properties  of  the  bulk 
of  such  transformed  specimens  but  would  be  strongly  biased  by  the 
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properties  of  the  g-SiC  "skin".  The  present  study  was  therefore  | 

limited  to  the  latter  two  surfaces  preparations  where  specimens  | 

with  uniform  microstructures  could  be  prepared  by  grinding  from  j 

3 

i 

larger  sintered  bodies.  i 

1 

A series  of  100  specimens  were  prepared  in  a manner  similar  i 

to  that  described  in  Section  II.  The  bars  were  sintered  in  five 
groups  of  20  at  temperatures  of  2055,  2075  (2  groups),  2095  and 
2115°C  for  30  minutes  to  produce  an  acceptable  range  of  transformation 
with  a minimum  variation  in  density  (95-96%  of  theoretical) . The 
as-sintered  dimensions  were  approximately  0.165  x 0.165  x 1.65  in. 

In  order  to  avoid  the  gradient  in  microstructure  described  above,  each 
specimen  was  diamond  ground  to  a cross  section  of  0.100  x 0.100  in. 

The  relative  amounts  of  3 and  a polytypes  were  determined  by 

4 

X-ray  diffraction  techniques.  The  results  are  included  in  Figure  3 
as  percent  of  transformation  to  ot-SiC  polytypes  versus  the  sintering 
temperature.  The  amount  of  transformation  varied  from  approximately 
5 percent  at  the  lowest  sintering  temperature  to  approximately  95  per- 
cent at  the  highest  temperature.  Two  of  the  five  groups  were  sintered 
at  2075®C  (separate  firings)  but  resulted  in  significantly  different 
amounts  of  transformation  as  seen  in  Figure  3. 

The  microstructures  of  the  five  firings  are  included  in  Figures  4a- 
e as  relatively  low  magnification  optical  micrographs  of  polished  and 
etched  sections.  The  polished  sections  were  etched  with  boiling  Murakami's 

reagent  (60  g KOH  + 60  g KjFe(CN)g  + 120  ml  H2O)  for  five  minutes.  As 

3 1 

described  elsewhere,  this  etch  reveals  primarily  a-g  and  a-a  grain  ! 

boundaries  and  only  reveals  g-g  grain  boundaries  as  very  low  contrast 

J 
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PERCENT  OF  TRANSFORMATION 


( 


( 


SINTERING  TEMPERATURE  ('O 

Figure  3.  Degree  of  p co  a Polytype  Transformation  Versus  Sintering 
Temperature  for  30  Minutes  at  Temperature 
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features  at  high  magnification  (molten  KCW  is  generally  more  useful 
for  revealing  6-g  grain  boundaries).  The  transformed  o-SiC  in 
these  microstructures  appears  in  two  distinct  morphologies  of 
"platelets"  and  "feathers",  both  of  which  have  been  described  in 
some  detail  in  References  1 and  3.  At  low  sintering  temperatures, 
platelets  appear  to  predominate  as  the  preferred  morphology  of  a 
growth,  while  feathers  predominate  at  higher  temperatures. 

Strength  testing  was  carried  out  at  room  temperature  in  three- 

point  bending  on  a one  inch  span  using  a sphere-loaded  testing 

2 

fixture  as  described  previously.  Of  each  group  of  20  specimens, 

10  were  tested  in  the  as-machined  state  and  10  were  refired  in 
vacuum  at  1960®C  for  40  minutes  before  testing.  As  discussed  in 
Section  III  this  refiring  treatment  is  generally  found  to  decrease 
the  severity  of  surface  flaws  such  as  those  produced  during  machining. 

The  relatively  low  refiring  temperature  of  1960®C  was  chosen  to 
minimize  the  risk  of  any  further  transformation  during  the  refiring 
treatment . 

The  results  of  bend  testing  are  included  on  Figure  5 as  fracture 
strength  versus  percent  of  transformation  to  a-SiC.  The  open  and  closed 
data  points  represent  the  average  strengths  of  as-machined  and  refired 
specimens  respectively  and  the  error  bars  include  the  95  percent  confidence 
limits  on  the  averages.  In  this  study,  no  significant  difference  in  strength 
could  be  resolved  between  the  specimens  tested  in  the  as-machined  condition 
and  those  tested  after  refiring.  This  lack  of  strength  increase  after 
refiring  suggests  that  the  strength  was  controlled  by  internal  flaws 
rather  than  machining  damage.  Based  on  the  conclusion  that  refiring  has 
no  significant  effect  on  the  strength  of  these  specimens,  a single  curve 
was  drawn  through  all  the  data  of  Figure  5.  The  strength  decreased  with 


FRACTURE  STRENGTH  (kpsi) 
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increasing  transformation  from  a maximum  of  approximately  85  kpsi 
(580  MPa)  at  100  percent  3 to  a minimum  of  approximately  60  kpsi 
(420  MPa)  at  100  percent  a-SiC. 


P 


It  should  be  emphasized  that  the  data  of  Figure  5 is  meaningful 
only  for  SiC  which  transforms  into  microstructures  similar  to  those 
shown  in  the  micrographs  of  Figure  4.  Transformation  into  other 
microstructures  such  as  small  platelets  or  uniform,  equiaxed  grains 
of  a-SiC  would  undoubtedly  result  in  a different  behavior  of  strength 
versus  transformation. 


a.  2055°C 

Figure  4.  Polished  and  Etched  Sections  Illustrating 
the  Effect  of  Sintering  Temperature 
on  Microstructure.  lOOX.  Etched 
with  Murakami  reagent. 
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b.  2075° 

Figure  4.  (Cont'd) 


d,  2095®C 
Figure  4.  (Cont'd) 


e.  2115°C 
Figure  4.  (Cont'd) 


V.  EFFECT  OF  OXIDATION  ON  FRACTURE  STRENGTH 


I 


The  successful  application  of  non-oxide  ceramics  as  high 
temperature  structural  materials  will  require  adequate  resistance 
to  oxidation  and  corrosion  under  operating  conditions.  Both  SiC 
and  SijN^  are  thermodynamically  unstable  in  oxidizing  environments, 
but  the  rates  of  oxidation  are  generally  controlled  by  the  mobility 
of  oxygen  through  the  oxide  layer.  In  cases  of  theoretically  dense, 
chemically  pure  SiC's  and  Si^N^'s  (such  as  those  made  by  CVD  processes) 
the  rates  of  oxidation  are  remarkably  slow,  even  at  temperatures  of 
1550  to  1600®C,^’^  thus  illustrating  the  excellent  potential  of 
pure  Si02  as  a barrier  to  oxygen  diffusion.  Spalling,  cracking  or 
chemical  contamination  of  the  Si02  layer  will  obviously  increase 
the  mobility  of  oxygen  through  the  oxide  layer  and  are  therefore 
undesirable. 

The  oxidation  behavior  of  boron-doped,  sintered  SiC  has  been 
reported  at  rates  comparable  to  that  of  CVD  SiC.^  Such  behavior  can 
probably  be  attributed  to  the  relatively  high  purity  of  the  oxide  layer, 
especially  if  the  1/2  w/o  boron  addition  is  continually  removed  from  the 
Si02  as  volatile  820^.  In  light  of  the  excellent  oxidation  resistance 
of  sintered  SiC,  a relatively  low  priority  was  assigned  to  the  task 
of  characterizing  the  fracture  strength  as  a function  of  exposure  to 
oxidizing  environments.  This  priority  was  raised,  however,  by  recent 
reports  in  the  literature  of  unexpectedly  large  strength  losses  after 
oxidation  in  commercially  available  hot-pressed  SiC's  and  Si^N^'s.^'^ 

An  an  example.  Ford  has  reported  in  their  ARPA  contract^  that  the 
average  bend  strength  of  Norton's  hot-pressed  Si^N^,  NC  132,  decreased 
from  110  to  35  kpsi  (760  to  240  MPa)  after  an  exposure  of  1000  hours 


20 


at  only  1370*C. 


In  response  to  this  increased  concern,  a series  of  60  bend 
specimens  were  prepared  in  a manner  similar  to  that  described  in 
Section  II  to  study  the  response  of  sintered  SiC  to  static  oxidation 
in  air  at  1500“C  for  times  to  2000  hours.  The  specimens  were  sintered 
in  a carbon  resistance  furnace  at  approximately  2055‘’C  for  30  min. 
at  a furnace  pressure  of  0.4  torr,  resulting  in  final  densities  of 
approximately  94  percent  of  theoretical.  Prior  to  oxidation,  the 
specimens  were  ground  parallel  to  the  specimen  length  to  a cross 
sectional  dimension  of  0.150  x 0.150  in.  Twenty  of  the  bend  specimens 
were  reserved  as  controls,  while  the  remaining  40  specimens  were  oxidized 
in  a MoSi2  resistance  heated  furnace.  The  specimens  were  held  vertically 
by  placing  them  in  shallow  holes  drilled  in  a hot-pressed  billet  of 
ZrSiO^.  Zircon  was  chosen  for  the  specimen  holder  because  of  its  com- 
patibility with  SiC  and  Si02  in  oxidizing  environments.  The  furnace 
was  initially  loaded  with  all  40  specimens,  heated  to  1500®C  for  10 
hours  and  cooled  to  room  temperature,  after  which  10  specimens  were 
removed.  In  a similar  fashion,  the  remaining  30  specimens  were  removed 
in  3 groups  of  10  after  integrated  exposures  of  200,  1000,  and  2000  hours 
at  1500’’C. 

Four  of  the  ten  specimens  exposed  for  the  full  2000  hours  were 
weighed  initially  and  each  time  the  furnace  was  cooled,  thus  providing 
data  on  weight  gain  due  to  oxidation.  These  results  are  included  as 
the  open  circles  on  Figure  6 where  weight  gain  is  plotted  versus  the 
square  root  of  time.  The  circles  represent  the  averages  of  the  four 
weight  gain  measurements  and  the  error  bars  include  the  95  percent 
confidence  intervals  on  the  averages.  The  data  points  fall  quite 
well  on  a straight  line  with  a slope  of  0.054  mg/cm  /hr  , thus 
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suggesting  a parabolic  rate  of  oxidation.  Parabolic  kinetics  are 
conunonly  observed  when  the  rate  of  oxidation  is  controlled  by  the 
diffusion  of  oxygen  through  a coherent,  non-spalling  oxide  layer.  The 
appearance  of  the  oxide  on  these  specimens,  however,  suggested  that  some 
spalling  began  between  200  and  1000  hours  and  continued  through  the  end 
of  the  experiment.  In  order  to  estimate  the  total  amount  of  spalling 
after  2000  hours,  the  four  specimens  that  were  weighed  through  the 
experiment  were  placed  in  hydrofluoric  acid  to  remove  all  of  the  Si02 
and  were  weighed  once  again.  The  difference  between  this  weight  and  the 
initial  (pre-oxidation)  weight  of  each  specimen  represented  the  amount 
of  SiC  removed  by  oxidation.  Since  each  mole  of  oxidized  SiC  formed  one 
mole  of  Si02,  one  can  easily  predict  the  total  weight  gain  including  that 
Si02  which  spalled  from  the  surface.  The  single  square  data  point  on 
Figure  6 represents  this  calculated  total  weight  gain  due  to  oxidation 
after  2000  hours  at  1500®C.  The  difference  between  the  circle  and  the 
rquare  data  points  at  2000  hours  is  the  best  estimate  of  the  weight  of 
Si02  lost  due  to  spalling.  The  non-linear  dashed  curve  included  on 
Figure  6 represents  an  estimate  of  the  total  weight  gain  versus  time 
corrected  for  errors  due  to  spalling.  It  was  drawn  tangent  to  the 
parabolic  oxidation  line  at  times  less  than  200  hours  where  no  spalling 
was  observed  and  was  connected  to  the  data  point  at  2000  hours  with  a 
smooth  relationship. 

The  fracture  strengths  of  all  60  specimens  were  measured  in  three- 
point  bending  at  room  temperature  using  a sphere-loaded  test  fixture  with 
a 1.5  in. span.  The  results  are  included  on  Figure  7 as  fracture  strength 
versus  the  square  root  of  exposure  time  at  1500®C.  (This  was  a convenient 
time  axis  to  plot  the  range  of  exposure  times  including  zero.)  Again, 
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the  data  points  represent  the  average  strengths  and  the  error  bars 
include  the  95  percent  confidence  intervals  on  the  averages.  All 
of  the  average  strengths  were  found  to  fall  between  60  and  80  kpsi 
(440  and  530  MPa).  The  dashed  line  drawn  through  these  data  requires 
some  explanation  because  there  is  more  information  available  than 
that  which  is  shown  on  Figure  7.  First,  considering  only  the  20 
control  specimens  and  the  10  specimens  exposed  at  1500®C  for  10 
hours,  the  distribution  of  the  difference  between  means  using 
Student's  distribution  reveals  that  there  is  approximately  a 90 
percent  chance  that  the  10  hour  specimens  come  from  a population 
with  a higher  mean  strength  than  the  as-machined  control  specimens. 

This  probability  increases  even  further  as  one  groups  the  10  and 
200  hour  specimens  together,  etc.  It  is  speculated  that  the  cause 
of  this  probable  strength  increase  may  be  another  "crack  blunting" 
mechanism  where  surface  defects  are  reduced  in  severity  by  material 
removal  at  the  crack  tips,  in  this  case,  by  oxidation.  Similar 

g 

observations  have  been  reported  by  Lange  for  hot -pressed  and 
reaction-bonded  SiC's.  Since  no  data  were  taken  between  0 and  10 
hours  of  exposure,  the  details  of  the  strength  versus  time  relationship 
are  unknown,  therefore  a dashed  straight  line  was  drawn  schematically 
between  the  two  points. 

The  average  strengths  of  Figure  7 between  10  and  1000  hours  appear 
to  plateau  at  approximately  76  kpsi  (525  MPa)  but  drop  to  65  kpsi  (450  MPa) 
by  2000  hours  of  exposure.  Fractographic  information  suggests  that  this 
drop  in  strength  is  indeed  real  but  is  probably  not  as  abrupt  as  suggested 
by  the  average  strengths.  The  fracture  origins  of  all  60  specimens  were 
located  and  cataloged.  In  the  case  of  those  oxidized  for  1000  hours, 

3 or  possibly  4 of  the  10  specimens  were  found  to  fail  from  a type 
of  surface  "pit"  which  was  not  evident  in  specimens  oxidized  for  shorter 
times.  In  the  specimens  oxidized  for  2000  hours,  these  pits  were  larger 


in  size  and  accounted  for  at  least  6 of  the  10  failures.  Removal 
of  the  Si02  from  several  specimen  fragments  revealed  that  the  surfaces 
contained  many  localized  areas  which  were  undergoing  greatly  accelerated 
oxidation.  (Due  to  orientation  effects,  etc.,  the  areas  of  localized 
attack  appear  to  be  related  to  the  cold  pressing  step  of  the  specimen 
fabrication) . Figures  8a-e  are  relatively  low  magnification  optical 
micorgraphs  (using  oblique  lighting)  of  the  oxidized  surfaces  with 
all  SiO^  removed.  The  five  micrographs  correspond  to  0 (as-machined), 

10,  200,  1000  and  2000  hours  of  exposure  at  1500°C  respectively.  The 
initial  stages  of  localized  attack  are  evident  even  at  10  hours  of 
oxidation  (Figure  8b) , but  the  resulting  pits  are  apparently  not 
sufficiently  severe  defects  to  cause  failure.  After  1000  and  2000 
hours  of  oxidation  (Figures  8d  and  e) , the  larger  pits  measure  50  to 
lOOvim  in  size  and  are  sufficiently  severe  to  decrease  the  average  frac- 
ture strength.  The  dashed  curve  on  Figure  7 beyond  10  hours  of  exposure 
was  drawn  schematically  as  a straight  line  and  was  included  primarily 
to  show  that  the  strength  decreased  with  increasing  exposure  time  in 
some  fashion  that  is  not  yet  satisfactorily  defined. 

In  summary,  the  rate  of  static  oxidation  of  sintered  SiC  in  air 
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at  1500®C  is  nearly  parabolic  with  a rate  constant  of  0.054  mg/ cm  /hr  ' . 
The  first  few  hours  of  oxidation  appear  to  increase  the  average  strength 
approximately  10  to  15  percent,  possibly  through  a crack  blunting  me- 
chanism. An  accelerated  mode  of  localized  oxidation  was  found  to  cause 
surface  pitting  and  an  associated  decrease  in  fracture  strength  after 
exposures  approaching  2000  hours.  Nevertheless,  the  average  strength 
of  specimens  tested  in  bending  was  no  lower  after  2000  hours  of  exposure 
at  1500®C  than  the  strength  of  as-machined  specimens. 


Surfaces  After  Oxidation  at  1500°C  and 
Removal  of  Si02  in  HF 


b.  10  hrs 


Figure  8.  (Cont'd) 


VI . FRACTURE  STATISTICS 

The  ability  to  effectively  design  highly  stressed  con5>onents 
of  any  structural  ceramic  will  require  an  accurate  description  of 
the  statistical  nature  of  fracture.  Unlike  most  metals,  ceramics 
will  require  this  additional  design  criteria  to  account  for  the 
large  variability  in  strength  which  is  commonly  observed  in  brittle 
materials  and  to  account  for  the  size  dependence  of  strength  which  is 
a natural  consequence  of  the  variability.  A proper  description  must 
relate  the  stress  at  which  failure  occurs  to  several  dependent  variables 
including: 

• the  cumulative  probability  of  failure — higher  stresses 
can  be  supported  if  higher  probabilities  of  failure  are 
acceptable; 

• the  specimen  size — larger  specimens  have  a greater 
probability  of  containing  a large  flaw  and  therefore 
can  support  less  stress;  and 

• the  stress  distribution — a body  under  uniform  tensile 
stress  has  a higher  probability  of  failing  at  its  most 
deleterious  flaw  than  a similar  body  with  only  a small 
fraction  of  its  volume  under  the  maximum  stress . 

The  probabilistic  nature  of  fracture  in  brittle  materials  originates 
from  variability  in  the  severity  of  defects.  In  turn,  the  statistical 
functions  which  are  used  to  describe  the  macroscopic  effects  of  such 
variability  can  be  reduced  to  descriptions  of  the  flaw  size  distribution. 
When  viewed  in  this  manner,  it  is  evident  that  a proper  description  of 
the  statistical  aspects  of  fracture  in  a given  material  is  extremely 
sensitive  to  any  factors  in  the  material  fabrication  process  which  may 
affect  that  flaw  size  distribution.  Such  factors  certainly  include  changes 
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in  powder  processing  and  material  fabrication  techniques.  As 

will  be  described  in  more  detail  in  Section  IX,  the  large  internal 

voids  in  sintered  SiC  (originally  discussed  in  last  year's  final 
2 

report  ) are  probably  caused  by  weak  agglomerates  formed  during 
the  powder  processing  steps.  Work  is  under  way  in  a concurrant 
program  to  identify  modified  powder  processing  and  material  fabrication 
techniques  to  avoid  such  defects.  In  light  of  this  potential  change 
in  both  the  powder  processing  technique  and  the  resulting  flaw  size 
distribution,  the  decision  was  made  to  delay  the  massive  testing  program 
necessary  to  satisfactorily  characterize  the  fracture  statistics  of 
the  present  form  of  sintered  SiC.  Nevertheless,  efforts  in  the  area 
of  fracture  statistics  have  continued  during  the  last  year  but  were 
directed  toward  the  development  and  verification  of  both  e^erimental 
and  mathematical  tools  to  be  used  in  the  future.  Some  initial  spin 
testing  has  been  carried  out  on  sintered  SiC  and  the  statistical 
parameters  of  all  the  various  groups  of  bend  specimens  have  been 
calculated,  but  the  two  areas  of  fracture  statistics  to  be  discussed 
in  detail  here  are: 

• A means  of  quantitatively  determining  the  statistical 
parameters  of  fracture  from  the  positions  of  fracture 
origins  was  derived  and  demonstrated  for  two-parameter 
Weibull  statistics. 

• TVo  sets  of  literature  data  were  reanalyzed  to  quanti- 
tatively illustrate  examples  of  multiple  flaw  distri- 
butions. 

The  most  useful  mathematical  descriptions  of  fracture  statistics 
presently  available  are  based  on  weakest  link  models  combined  with 


I • empirical  distribution  functions.  One  very  popular  example  is  Weibull 
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^ fracture  statistics  which  is  generally  reported  to  describe  experimental 

data  reasonably  well.  It  should  be  emphasized,  however,  that  the 

a/ 

Weibull  distribution  function  is  not  derived  completely  from  first 
principles  but,  rather,  is  an  empirical  description  of  the  probability 
of  failure.  In  fitting  experimental  data  to  such  a mathematical 
j function,  a great  deal  of  caution  must  be  exercised  in  extrapolating 

the  function  beyond  the  data  base.  In  spite  of  the  weaknesses  caused 
by  its  empirical  nature,  the  relatively  simple  distribution  function 
of  Weibull  fracture  statistics  provides  a useful  model  to  illustrate 
mathematical  manipulations  that  may  require  numerical  methods  for  a 
more  complex  distribution  function.  For  this  reason,  the  Weibull 
distribution  was  used  in  both  of  the  discussions  to  follow.  Since 
both  sections  will  refer  to  specifics  of  the  Weibull  distribution 
function,  it  will  be  useful  to  begin  by  reviewing  the  form  of  the 
I distribution  function. 

The  Weibull  distribution  function  relates  the  cumulative  pro- 
bability of  failure  (P)  to  the  volume  (V)  under  tensile  stress  (o) 
with  two  parameters- -the  Weibull  modulus  (m)  and  a normalizing  constant 
(Oq)  as  follows: 

P 

This  form  of  the  function  has  two  adjustable  parameters  and  integrates 

over  the  volume  of  the  specimen.  Thus  it  is  generally  referred  to  as  ; 

two-parameter,  volume  Weibull  statistics.  Other  variations  of  the 

Weibull  distribution  function  can  include  a third  ajustable  parameter 

and  may  integrate  over  the  surface  of  the  body  rather  than  the  volume  ; 

hmJ ^ 
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(if  all  flaws  are  expected  to  be  surface  flaws).  The  Weibull 
modulus,  m,  is  a measure  of  the  scatter  or  dispersion  of  the 
distribution  (similar  to  the  standard  deviation  of  a normal  dis- 
tribution) where  a small  m represents  a large  degree  of  scatter. 

Most  ceramics  are  reported  to  have  m's  in  the  range  of  four  to  ten, 
whereas  metals  are  estimated  to  have  m's  in  the  range  of  50  to  100. 

The  normalizing  constant  a^,  is  an  indirect  measure  of  the  average 
strength.  For  instance,  if  the  Weibull  modulus,  the  specimen  size 
and  the  loading  configuration  are  held  constant,  the  average  strength 
is  directly  proportional  to  o^. 

For  boi\es  under  uniform  tension,  the  integration  of  Equation 
1 is  trivial,  but  for  the  general  case  where  the  stress  in  a body 
is  a function  of  position,  the  integration  must  be  carried  out 
either  analytically  or  numerically.  As  an  example,  the  stress  in  a 
rectangular  beam  in  three-point  bending  can  be  described  analytically 
as: 

a = 6 P X y (2) 

w h^ 

where  P is  the  center  point  load,  x is  the  distance  from  the  nearest 
outer  load  point,  y is  the  distance  from  the  neutral  axis,  w is  the 
specimen  width  and  h is  the  specimen  height.  Substituting  this  into 
Equation  1,  expanding  the  volume  integral  into  a triple  integral  and 
rearranging  results  in: 


P = 1-exp 
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where  js  the  maximum  tensile  stress  in  the  beam.  Integration  of 
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this  expression  yields: 


Li 


P = l-exp  I j 

I 2(m+l)^ 


''  (-v)l 


P = l-exp  I - 


where  K is  a "load  factor" 


V - ^ 

2 

2(m+l) 

The  load  factor  for  the  more  general  case  of  four-point  bending  with 
volume  flaws  can  be  derived  from  Equation  1 in  a similar  manner  and 
results  in: 

(J)m  + 1 

Kv  = ^ ; (7) 

2Cm+l)2 

where  a is  the  inner  span  and  b is  the  outer  span.  Equation  7 cor- 
rectly reduces  to  Equation  6 for  three  point  bending  where  the  inner 
span,  a,  is  zero.  In  those  cases  where  failure  is  controlled  by 
surface  flaws.  Equation  1 should  be  integrated  over  all  elements  of 
surface  rather  than  volume.  The  resulting  load  factor  for  four-point 
bending  can  be  related  to  Equation  7 as  follows: 

•'a  = 

where  w and  h are  the  specimen  width  and  height  respectively. 

The  product  of  Ky  times  V is  often  referred  to  as  the  "effective 
volume"  under  tension  and,  as  the  term  implies,  can  be  pictured  as 
the  volume  of  material  which  is  effectively  being  tested  under  uniform 
tension.  Accordingly,  times  surface  area  is  the  effective  area  under 
uniform  tension.  Equation  1 can  be  integrated  either  analytically  or 
numerically  to  the  form  of  Equation  S for  a stressed  body  of  any 


geometry  and  loading  configuration.  Except  for  the  case  of  uniform 

tension,  the  resulting  load  factors  will  always  be  numerically  less 

than  one  and  will  always  be  a function  of  m. 

The  Weibull  modulus  and  a can  be,  estimated  from  fracture  data 

o 

by  a number  of  methods,  but  all  of  the  common  ones  can  be  classified 
into  one  of  two  categories.  The  most  popular  methods  require  fracture 
strength  measurements  on  a number  of  specimens  of  identical  geometry. 

The  fracture  stengths  are  ordered  from  weakest  to  strongest,  and  each 
is  assigned  a probability  of  failure  based  on  its  ranking,  n,  as  follows; 


(9) 


where  N is  the  total  number  of  data  points.  A regression  analysis  can 


then  be  performed  using  Equation  5 at  constant  KV  to  find  the  best  values 
for  m and  to  describe  all  N data  points  (each  data  point  consists  of  a 
fracture  strength  and  an  associated  probability  of  failure) . Variations 
of  this  method  arise  from  different  possible  weighting  factors  of  the 
data  in  the  regression  analysis. 


One  of  these  methods  is  very  convenient  because  it  allows  a single 
linear  regression  analysis  of  the  data.  Equation  5 can  be  manipulated 
and  rearranged  into  the  following  form: 

In  In  (y^)  = m Ina^^^  + ln(  *'^— ) (10) 

This  is  an  equation  of  a straight  line  and  is  therefore  easily  subjected 
to  linear  regression.  Equation  10  also  allows  a convenient  graphical 
representation  of  the  data  as  shown  schematically  in  Figure  9 by  plotting 
In  In  ("Y^p  ) vs  In  The  slope  of  the  least  squares  line  through  the 

data  is  m and  the  position  of  the  line  defines  It  should  be  noted 

that  the  maximum  stress  of  each  specimen  is  used  in  all  of  this  discussion 
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despite  the  fact  that  fracture  in  three-point  bend  specimens,  for 
example,  virtually  never  originates  exactly  at  the  line  of  maximum 
tensile  stress. 


The  other  common  method  of  estimating  the  Weibull  parameters 
requires  the  measurement  of  fracture  strength  at  a variety  of 
effective  volumes.  Rearranging  Equation  5 again  yields: 

In  cr„  ^ In  KV  + In  o + — — in  In  — (11) 

Mm  o m 1-P 

At  a constant  probability  of  failure,  this  defines  a linear  relation- 
ship between  and  In  KV.  Graphical  representation  of  this  equation 
is  shown  schematically  in  Figure  10  as  In  versus  In  KV  with  the 
slope  of  the  line  equal  to  — The  data  to  be  plotted  on  such  a 
graph  should  be  at  constant  probability  of  failure,  thus  one  can 
plot  average  strengths  or,  somewhat  more  correctly,  median  strengths. 
Again,  variations  of  this  method  arise  from  different  types  of 
regression  analysis  on  the  data  points. 

I Graphical  representations  of  data  such  as  Figures  9 and  10 

also  allow  judgements  concerning  the  quality  of  the  distribution 
fimction  in  describing  the  data.  If  the  data  points  are  clearly 
non-linear  on  such  plots,  then  the  distribution  fimction  does  not 
adequately  describe  the  behavior,  regardless  of  the  choice  of 
statistical  parameters  m and  a^. 

VI-A.  Derivation  of  Statistical  Parameters  of  Fracture  from  the 
Locations  of  Fracture  Origins 

The  previous  review  of  fracture  statistics  has  been  included 
as  background  on  which  to  base  the  derivation  of  yet  another  in- 
I dependent  means  of  estimating  the  statistical  parameters  and  judging 
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in  cr 

Figure  9.  Schematic  Drawing  Illustrating  a Means  of 


Estimating  the  WeibuU  Modulus  from  Prob- 
ability of  Failure  and  Fracture  Strength 


L06  EFFECTIVE  VOLUME  (kV) 


Figure  10.  Schematic  Drawing  Illustrating  a Means  of  Estimating  the 
WeibuU  Modulus  from  Fracture  Strength  and  Effective 
Volume 
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the  validity  of  a distribution  function.  This  method  uses  additional 

! 

information  from  fracture  specimens  tdiich  is  easily  measured  but  not 
normally  compiled.  The  method  is  restricted  to  test  specimens  which 
contain  stress  gradients  (such  as  bend  specimens)  and  the  additional  | 

measurements  necessary  for  the  analysis  are  the  positions  of  the 
fracture  origins . 

The  distribution  of  positions  of  fracture  origins  is  strongly 
influenced  by  the  variability  in  strength  as  can  be  demonstrated  by 
first  considering  the  limiting  case  of  a material  with  no  variability 
in  strength  (standard  deviation  = 0;  m = <») . In  such  a hypothetical 
material,  failure  will  always  occur  at  exactly  the  same  stress,  and 
no  element  of  volume  can  survive  a higher  stress.  When  tested  in 
I three-point  bending  the  fracture  origins  will  always  be  located  exactly 

on  the  line  of  maximum  tensile  stress.  If  one  now  considers  a material 
with  some  variability  in  strength,  the  most  probable  position  for  failure 
I • in  bending  remains  on  the  line  of  maximum  tensile  stress,  but  not  all 

failures  will  initiate  there.  Due  to  the  var*-=': ility  in  strength  from 
location  to  location,  failure  will  often  occur  at  positions  where  the 
stress  is  somewhat  lower  than  the  maximum.  Accordingly,  the  fraction 
of  origins  located  away  from  the  line  of  maximum  tensile  stress  will 
increase  with  increasing  variability  in  fracture  strength  (or  decreasing 
m) , thus  providing  a means  of  estimating  the  variability  in  strength  by 
measuring  only  the  positions  of  fracture  origins. 

Two-parameter  volume  Weibull  statistics  and  three-point  bend 
specimens  have  been  chosen  to  illustrate  the  mathematics  of  this 
method,  but  the  underlying  logic  should  be  valid  for  any  distribution 
function  of  interest  and  any  loading  configuration  except  uniform 
tension.  The  position  of  a fracture  origin  in  a three-point  bend  ! 
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specimen  can  be  related  to  the  ratio  of  the  stress  at  the  fracture 
origin  to  the  maximum  stress  in  the  beam  without  any  knowledge  of  the 
actual  fracture  stress,  load,  etc.  by  the  following  relationship: 


where  x is  the  distance  from  the  fracture  origin  to  the  nearest  outer 


L.. 


load  point,  y is  the  distance  from  the  origin  to  the  neutral  axis,  h 
is  the  specimen  thickness  and  L is  the  specimen  length.  This  ratio 
of  stress  to  maximum  stress  will  be  referred  to  as  the  relative  stress, 
R.  Contours  of  constant  relative  stress  can  be  calculated  for  a beam 


by  rearranging  Equation  12  to  yield: 


R hi 
X 4 


For  a constant  value  of  relative  stress,  equation  describes 

a family  of  hyperbolae  with  a common  origin  on  the  neutral  axis  at  an 
outer  load  point.  Figure  11  is  the  side  view  of  a bend  specimen  showing 
several  contours  of  constant  relative  stress. 

The  total  probability  of  failure  of  a three-point  bend  specimen 
was  derived  earlier  by  integrating  Equation  3 with  appropriate  limits 
of  integration  to  include  the  entire  volume  of  the  specimen  under 
tensile  stress.  In  deriving  m from  the  positions  of  fracture  origins 
in  a material  with  volume  flaws,  a similar  integration  must  be  carried 
out  except  that  the  limits  of  integration  are  chosen  to  include  only 
that  portion  of  the  specimen  at  a relative  stress  greater  than  a given 
limit.  For  instance,  if  the  chosen  relative  stress  were  0.8,  then  the 
result  would  be  the  probability  that  fracture  would  originate  in  the 
volume  of  material  contained  within  the  0.8  relative  stress  contour  of 
Figure  11.  To  distinguish  it  from  the  total  probability  of  failure, 

P or  P.J,,  the  result  of  this  integration  will  be  referred  to  as  Pj^ 
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for  probability  of  failure  at  a relative  stress  greater  than  R.  The 
form  of  this  equation  for  three-point  bending  with  the  proper  limits 
of  integration  is: 

\m  u k/9  I/")  m 

dx  dy  dz 


1-exp 


[■‘  ^ 4-  / /■ . 


(xy)" 
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which  can  be  integrated  to: 

= 1-exp 

where : 


h (v)-] 


c = 1-r'"'"^  (1-ln  R®*^) 


(14) 


(15) 


(16) 


For  the  case  of  R=0,  the  integration  is  carried  out  over  the  entire  volume 
of  the  specimen,  therefore  Pj^  should  equal  P^.  Substituting  R=0  into 
Equation  16  results  in  C=l,  thus  Equation  IS  correctly  reduces  to 
Equation  5.  The  other  simple  limiting  case  is  R=l.  Since  there  is 
no  volume  at  a relative  stress  greater  than  one,  the  probability  of  failure 
in  that  volume  must  be  zero.  Substituting  R=1  into  Equation  16  results 


in  C*0,  thus  Equation  15  correctly  predicts  Pj^=0. 

More  useful  than  Pj^  is  the  ratio  of  Pj^  to  P^.  For  a given  value  of 
relative  stress,  this  is  the  probability  of  failure  occurring  in  positions 
of  the  specimen  within  that  stress  contour  divided  by  the  total  probability 
of  failure  of  the  entire  specimen.  This  ratio  therefore  corresponds  to 
the  fraction  of  fracture  origins  which  are  located  at  a relative  stress 
greater  than  R.  r I Im1 
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Figure  11.  Contours  of  Constant  Relative  Tensile  Stress  in  a Three -Point 
Bend  Specimen 


RELATIVE  STRESS  AT  THE  FRACTURE  ORIGIN  » R 


Figure  12.  Schematic  Drawing  Illustrating  a Means 
of  Estimating  the  Weibull  Modulus  from 

Pr/Pt  ^ 


If  this  ratio  is  independent  of  the  total  probability  of 


failure,  P^,  then  the  Weibull  modulus,  m,  can  be  predicted  solely 
from  the  positions  of  fracture  origins  with  absolutely  no  knowledge 
of  the  magnitude  or  variability  of  the  actual  fracture  strengths.  A 
graphical  approach  for  such  a situation  is  shown  in  Figure  12  where 
is  plotted  schematically  as  a function  of  R for  several  values 
of  m.  Experimentally  the  position  of  the  fracture  origin  in  each  of 
several  bend  specimens  would  be  measured  and  the  values  of  relative 
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stress  calculated  by  Equation  12.  These  relative  stresses  would  be 
ranked  in  order  of  decreasing  R from  1 to  N and  the  values  of  Pj^/Pj 
estimated  by: 

VPt=  nTT 

The  resulting  N data  points  would  then  be  plotted  on  Figure  12  (as 
shown  for  a hypothetical  set  of  nine  data  points)  and  the  Weibull 
modulus  estimated  from  the  positions  of  the  data  relative  to  the 
plotted  curves.  More  accurate  estimations  could  be  calculated,  of 
course,  by  iteratively  searching  for  the  value  of  m with  the  least 
scatter  of  data  about  the  predicted  line. 

Unfortunately,  the  relatively  simple  procedure  described  in 
the  previous  paragraph  cannot  be  used  because  the  ratio  of  Pn/P- 
is  not  independent  of  the  total  probability  of  failure.  Stated  another 
way,  the  ratio  is  a function  of  the  maximum  stress  at  the  time  of 
failure,  o , as  seen  on  the  right  side  of  Equation  17.  This  de- 
pendence  of  Pj^/Py  on  the  total  probability  of  failure  is  illustrated 
quantitatively  on  Figure  13  where  P„/P  is  plotted  as  a function  of 
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R (same  axes  as  Figure  12)  for  several  values  of  P.j,  at  a constant  Weibull 
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RELATIVE  STRESS  AT  THE  FRACTURE  ORIGIN  s R 

Figure  13.  Pr/Pt  Versus  R at  Constant  m Illustrating  the  Dependence 
of  the  Position  of  Fracture  Origins  on  the  Total  Probability 
of  Failure,  P-j- 
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modulus  of  7.  At  low  total  probabilities  of  failure,  the  fionction 


is  insensitive  to  P^,  but  as  the  total  probability  of  failure 
approaches  unity  there  is  a distinct  tendency  for  failure  to  occur 
closer  to  the  line  of  maximum  stress.  In  the  limit  of  = 1,  all 
failures  should  occur  exactly  on  the  line  of  maximum  tensile  stress. 
Stated  in  other  words,  the  probability  of  a specimen  failing  at  a 
given  position  is  not  only  a function  of  the  Weibull  modulus,  but 
also  depends  on  whether  that  particular  specimen  is  the  strongest 
or  weakest  specimen  of  the  batch. 

At  this  point,  there  are  four  interrelated  variables  which  are 
of  interest:  m,  R,  and  P^.  All  four  variables  cannot 

conveniently  be  illustrated  at  the  same  time  in  two  dimensions, 
therefore  Figures  14  and  15  both  plot  relative  stress  at  the  fracture 
origin  versus  total  probability  of  failure.  Figure  14  includes  a 
family  of  curves  where  Pj^/Pj  was  varied  and  the  Weibull  modulus 
held  constant  at  m = 7.  Figure  15  includes  a similar  family  of 
curves  where  m was  varied  and  Pj^/P^.  was  held  constant  at  0.5. 

The  simplified  method  discussed  above  (Figure  12)  required 
estimates  of  R and  Pj^/P.p  for  each  specimen,  after  which  m could  be 
predicted  graphically.  The  dependence  of  fracture  origin  position 
on  total  probability  of  failure  as  illustrated  in  Figures  13-15 
somewhat  complicates  this  method  of  predicting  m.  In  addition  to  R 
and  Pj^/Pj»  an  estimate  of  P^  is  needed  for  each  specimen.  Unfortunately, 
this  requires  some  information  concerning  the  actual  fracture 
strengths  which  was  not  needed  for  the  earlier  model.  The  ranking 
of  strengths  in  increasing  order  combined  with  Equation  9 is  sufficient 
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TOTAL  PROBABILITY  OF  FAILURE 

Figure  14.  R Versus  Prp  at  Constsuit  m for  Various 
Values  of  Pr/Pt 


total  probability  of  failure 

Figure  15.  R Versus  P-p  at  Constant  Pr/Pt  for  Various 
Values  of  m 


to  estimate  a total  probability  of  failure  for  each  specimen.  In 
order  to  estimate  the  Weibull  modulus,  the  resulting  data  are 
plotted  on  a three-dimensional  equivalent  of  Figure  12  (P.j,  is 
the  third  axis) • Since  three-dimensional  graphs  are  difficult  to 
construct,  the  most  feasible  approach  is  an  iterative  least  squares 
technique  with  the  aid  of  a computer. 

Experimental  verification  of  this  technique  required  a material 

with  only  a single  distribution  of  flaws  and  with  easily  identifiable 

2 

fracture  origins.  As  discussed  previously  , boron-doped,  sintered 
SiC  has  been  found  to  contain  several  different  flaw  distributions  and 
therefore  was  not  chosen.  Instead,  cylindrical  rods  of  pyrex  glass 
tested  in  three-point  bending  were  chosen  for  this  demonstration.  A 
series  of  45  specimens  of  5 mm  diameter  rod  were  cut  into  three  inch 
lengths.  To  insure  that  failure  occurred  from  only  one  flaw  distribution, 
the  specimens  were  uniformly  abraded  by  tumbling  for  two  hours  in  a 
slurry  containing  220  grit  SiC  abrasive  powder. 

The  previous  outline  for  estimating  m from  fracture  origins  was 
discussed  in  terms  of  a distribution  of  volume  flaws  in  a rectangular 
specimen.  The  analysis  therefore  had  to  be  revised  for  the  present 
case  of  a distribution  of  surface  flaws  in  a cylindrical  specimen.  The 
resulting  relationship  describing  Pj^  in  terms  of  R,  etc.  (the 
equivalent  to  Equation  14)  is  a double  integral  which  cannot  be 
completely  integrated. 


Pp  = 1-exp 


D is  the  specimen  diameter  and  0 is  the  angle  from  the  fracture  origin 
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to  the  position  of  maximum  tensile  stress.  The  remaining  integration 
can,  of  course,  be  performed  numerically. 

The  abraded  specimens  were  tested  at  room  temperature  in  three- 
point  bending  on  a two  inch  span  at  a loading  rate  of  approximately 
50  kpsi/min  (340  MPa/min)  and  resulted  in  an  average  strength  of  14.1 
kpsi  (97.2  MPa).  All  fracture  origins  were  located  on  the  specimen 
surface.  The  position  of  each  origin  was  defined  by  measuring  the 
distance  from  the  origin  to  the  nearest  outer  load  point  and  the 
angle  from  the  origin  to  the  point  of  maximum  tensile  stress.  These 
measurements  were  taken  on  each  half  of  the  specimen  and  averaged.  The 
relative  stresses  were  ranked  and  used  in  Equation  18  to  estimate  Pj^/P.j,, 
and  the  strengths  were  ranked  and  used  in  Equation  9 to  estimate  P.p. 

The  resulting  45  sets  of  R,  P.^,  and  Pj^/P-p  were  entered  into  a con^uter 
program  to  iteratively  find  the  most  suitable  value  of  m in  Equation  19 
to  describe  the  data.  The  result  was  a Weibull  modulus  of  12.9. 

For  the  sake  of  comparison,  the  Weibull  modulus  was  also  estimated 
in  a conventional  manner  from  the  fracture  strengths  using  the  linear 
regression  method  illustrated  earlier  in  Figure  9.  This  resulted  in  a 
Weibull  modulus  of  11.4  which  compares  quite  favorably  with  the  estimate 
of  12.9. 

In  summary,  a method  was  outlined  to  determine  the  Weibull  modulus 
of  a group  of  specimens  using  primarily  the  locations  of  fracture  origins 
in  specimens  with  stress  gradients.  The  method  does  not  employ  the  numerical 
values  of  fracture  stress  and  is  therefore  independent  of  those  methods  of 
estimating  m which  do  require  fracture  stresses.  The  exan^le  illustrated 
in  detail  was  two-parameter  volume  Weibull  statistics  on  a three-point  bend 
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specimen,  but  the  approach  which  was  outlined  can  be  applied  to 
different  specimen  geometries  and  distribution  functions.  This 
method  of  determining  the  Welbull  modulus  Is  not  meant  to  replace 
the  conventional  methods  reviewed  earlier,  but  rather  Is  meant  to 
supplement  them  by  more  effectively  using  all  of  the  Information  con- 
cerning the  variability  of  flaw  sizes  that  Is  available  In  a group  of 
test  specimens.  In  light  of  the  time,  effort  and  expense  associated 
with  the  fabrication  and  testing  of  fracture  specimens,  a means  of  nearly 
doubling  the  resulting  statistical  Information  per  specimen  Is  very  attractive. 
In  addition  to  providing  more  efficient  utilization  of  available  data,  the 
above  discussion  on  the  distribution  of  fracture  origin  positions  has 
yielded  useful  Information  on  two  related  subjects: 

• The  most  popular  view  of  fracture  statistics  attributes 
the  observed  effects  to  variations  in  flaw  size  from 
location  to  location  combined  with  the  weakest  link 
concept.  A less  popular  explanation  attributes  the 
effects  to  errors  in  the  testing  technique  including 
both  measurement  errors  and  errors  due  to  misalignments 
of  the  specimen  and  fixtures.  If  the  scatter  in  strengths 
is  primarily  due  to  testing  errors  instead  of  variations 
in  flaw  sizes,  then  there  should  be  little  or  no  dependence 
of  strength  on  specimen  size.  Also,  If  there  Is  Indeed 
no  variability  in  flaw  sizes  (the  true  m = <») , fracture 
should  always  originate  at  the  location  of  the  highest 
tensile  stress  and  the  only  scatter  in  fracture  origin 
positions  should  be  due  to  variability  In  the  nature 
of  the  parasitic  stresses  from  specimen  to  specimen.  It 
Is  very  unlikely  in  this  case  that  the  scatter  in  fracture 
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strengths  and  the  scatter  in  fracture  origin  positions 
would  coincidentally  lead  to  the  same  value  of  Weibull 
modulus.  Therefore  approximately  equivalent  values  of 
m as  calculated  by  the  two  techniques  (such  as  illustrated 
in  the  example  of  abraded  glass  rods)  helps  to  support 
variability  in  flaw  sizes  as  the  source  of  the  observed 
scatter  in  fracture  strengths. 

• As  mentioned  earlier  in  this  section,  the  fracture  strength 

of  a specimen  is  correctly  quoted  as  the  maximum  tensile  stress 
in  that  specimen  even  though  failure  in  bending,  for  example, 
virtually  never  occurs  exactly  at  the  position  of  maximum  stress. 

A number  of  investigators  reporting  in  the  literature  have  been 
concerned  with  the  use  of  as  the  best  description  of  the  fracture 
strength  and  instead  have  measured  the  location  of  the  fracture 
origin,  calculated  the  stress  present  at  that  location  at  the 
time  of  failure  and  incorrectly  used  that  number  to  calculate 
the  average  strength  and  the  Weibull  modulus.  This  approach 
leads  to  significantly  low  estimates  of  both  the  average  strength 
and  the  Weibull  modulus,  but  both  can  be  quantitatively  related 
to  the  correct  values  through  the  information  on  Figures  IS-l.S. 

VI-B.  Multiple  Flaw  Distributions 

Careful  fractography  of  most  ceramic  materials  will  reveal  more 
than  one  type  of  fracture  initiating  flaw  within  a set  of  specimens. 

The  source  of  these  flaws  are  many  and  include  machining  damage,  in- 
clusions of  foreign  material,  large  voids  or  void  clusters  due  to 
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powder  agglomeration,  etc.  Each  type  of  flaw  is  present  in  the 

material  as  a distribution  of  flaw  sizes  and  densities  which,  in 

most  cases,  is  totally  independent  of  all  other  flaw  distributions. 

The  consequences  of  multiple  independent  flaw  distributions  was 

2 

discussed  in  some  detail  in  last  year’s  final  report,  but  no  data 
were  presented  to  substantiate  the  conclusions.  The  primary 
material  being  studied  in  this  contract  is  boron-doped  sintered  SiC 
and,  as  also  discussed  in  last  year's  final  report,  at  least  four 
different  flaw  distributions  are  active  in  three-point  bend  specimens 
of  this  material.  A very  large  number  of  specimens  would  be  necessary 
to  satisfactorily  characterize  each  flaw  distribution.  Such  information 
will  eventually  be  required  but  at  present  is  not  available  for 
sintered  SiC. 

The  effect  of  multiple  flaw  distributions  can  be  illustrated 
at  present,  however,  by  using  published  data  on  other  ceramic 
materials  which  are  reasonably  well  characterized  and  which  do  not 
exhibit  such  a large  number  of  active  flaw  distributions.  Two 
examples  of  such  materials  have  been  chosen  for  the  present  discussion. 
In  both  cases  the  original  authors  kindly  provided  more  detailed 
test  results  than  those  available  in  the  literature.  The  first 
example  is  based  on  a study  of  CVD  SiC  fibers  tested  in  uniaxial 
tension  by  Kotchick,  Hink  and  Tressler^^where  the  strength  of  a 
large  number  of  fibers  was  measured  as  a function  of  gage  length  at 
a constant  fiber  diameter.  Unfortunately,  due  to  fragmentation 
after  the  initial  specimen  failure,  it  was  not  feasible  to  locate 
and  classify  the  fracture  initiating  flaws.  The  second  example  to 
be  discussed  is  based  on  a study  of  a commercial  96%  ^^2^3 

11  12 

in  three  and  four-point  bending  by  Bansal,  Duckworth  and  Niesz,  ’ 
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where  the  strengths  of  a moderate  number  of  specimens  were  measured 
as  a function  of  specimen  size,  loading  configuration  and  environment. 

In  this  case  the  fracture  origins  of  all  specimens  were  identified, 
thus  helping  to  substantiate  the  existence  of  multiple  flaw  distributions. 
The  authors  of  both  examples  concluded  that  conventional  Neibull 
statistics  could  not  describe  the  observed  behavior.  The  present 
re-analysis  of  their  data  illustrates  that  two-parameter  Weibull 
statistics  describe  the  data  satisfactorily  (exceptionally  well  in 
the  case  of  Bansal's  data)  if  the  materials  are  assumed  to  each 
contain  two  independent  flaw  distributions. 

The  problem  of  estimating  the  applicability  of  a particular 

mathematical  ftinction  in  describing  the  probabilities  of  failure  of 

specimens  of  different  sizes,  geometries  and  loading  configurations 

quickly  reduces  to  a problem  of  graphically  visualizing  both  the 

data  and  the  predicted  behavior  simultaneously.  In  the  case  of  a 

set  of  specimens  of  similar  size,  geometry  and  loading  configuration, 

this  can  be  done  conveniently  by  plotting  the  data  as  In  In  yip-  vs. 

In  fracture  strength  as  shown  schematically  on  Figure  9.  If 

two-parameter  Weibull  statistics  are  correct,  the  data  should  fall 

on  a straight  line,  but  if  two  such  distributions  are  active,  the 

data  will  fall  on  intercepting  straight  lines  as  shown  schematically 

2 

in  Figures  27  and  28  of  last  year's  final  report.  Fracture  data 
derived  from  specimens  of  more  than  one  size  or  tested  in  more  than 
one  loading  configuration  cannot  all  be  simultaneously  plotted  in 
such  a manner,  therefore  one  must  manipulate  the  data  in  order  to 
satisfactorily  incorporate  all  of  it  onto  one  two-dimensional  plot. 
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Weibull  two -parameter  statistics  have  been  chosen  again  for 


1 


this  discussion  because  the  relatively  simple  function  allows 
mathematical  manipulations  which  are  not  possible  for  more  complex 
statistical  functions  and  because  it  very  nicely  describes  the  data 
of  Kotchick  et  al.  and  Bansal  et  al.  As  discussed  earlier.  Equation  5 
relates  the  maximum  applied  stress  and  the  specimen  volxnne  to  the 
probability  of  failure  for  a material  with  a single  distribution 
of  volume  flaws.  The  equivalent  relationship  to  Equation  5 for 
two  independent  volume  distributions  is: 


P = 


(20) 


where  the  subscripts  1 and  2 for  m's  and  o^'s  distinguish  the  parameters 
for  the  two  distributions.  Since  the  load  factors  are  normally  a 
function  of  m,  the  quantity  cannot  be  factored  out  of  the  sum  in 
the  exponential. 

In  the  case  of  fracture  data  from  multiple  specimen  sizes,  one 
must  choose  a single  specimen  size,  geometry  and  loading  configuration 
to  normalize  all  other  data.  Each  original  piece  of  fracture  data 
consists  of  a fracture  strength  and  an  estimated  probability  of 
failure.  Normalizing  a data  point  to  a different  specimen  size  can 
be  accomplished  in  two  ways;  The  strength  can  be  shifted  by  an 
appropriate  amount  with  no  change  in  probability  of  failure;  or  the 
probability  of  failure  can  be  shifted  with  no  change  in  the  strength. 

The  latter  of  these  two  methods  is  the  only  practical  choice  for 
multiple  flaw  distributions. 
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For  the  general  case  of  normalizing  data  from  one  arbitrary  specimen 
size,  geometry  and  loading  configuration  (designated  by  subscripts 
a)  to  another  (designated  by  subscripts  g) , one  can  write  two 
equations,  each  of  the  fonn  of  Equation  20,  as  follows: 


= 1 - exp 


Pg  = 1 - exp 


)1 


(21) 


['.  {■■■  • >»  -4>i 


(22) 


where  is  the  same  maximum  stress  in  both  specimens  and  where  both 
distributions  are  present  as  volume  flaws.  For  the  case  of  uniform 
uniaxial  tension  (Ky=l)  these  equations  can  be  combined,  rearranged  and 
simplified  to  yield: 


Pg  = 1 - (1 


(^) 

I V •' 


(23) 


It  is  interesting  to  note  that  the  shift  in  probability  of  failure 
with  a shift  in  volmne  (at  a constant  failure  stress)  is  totally 
independent  of  the  statistical  parameters,  m and  o^.  In  fact,  for 
the  case  of  constant  diameter  specimens  tested  in  uniaxial  tension 
(such  as  the  data  of  Kotchick  et  al.),  it  can  be  shown  that  the 
shift  in  probability  is  only  dependent  on  the  assumption  that 
weakest  link  statistics  are  applicable,  llte  amount  of  shift  is 
independent  of  the  specific  function  (two-parameter,  three-parameter 
Weibull,  etc.)  and  the  number  of  active  flaw  distributions  and  is 
also  correct  when  both  volume  and  area  distributions  are  active 
simultaneously. 
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A more  general  description  than  Equation  23  is  needed  for  the 
data  of  Bansal's  where  the  normalization  of  the  data  must  accommodate 
changes  in  specimen  size,  changes  in  loading  configuration  and 
active  distributions  of  both  area  and  volume  flaws.  The  equivalent 
relationship  to  Equation  23  for  this  general  case  is  more  complex: 


where. 
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In  this  case,  m and  0^  must  be  known  or  estimated  for  each  distribution 
before  any  of  the  data  can  be  shifted,  and  the  result  is  subject  to  any 
errors  in  the  parameters. 

10 

The  work  on  SiC  fibers  by  Kotchick,  Hink  and  Tressler  included  uni- 

* 

axial  testing  of  approximately  300  fibers  of  Avco  CVD  SiC  with 
seven  different  gage  lengths  from  0.25  to  20  in.  The  fibers  were 
approximately  0.004  in.  diameter  with  a carbon  core  of  approximately 
0.001  in.  diameter.  Although  SiC  fibers  were  also  tested  in  the 
abraded  condition,  the  300  fibers  of  interest  in  this  discussion 


were  tested  with  virgin  surfaces.  The  average  strengths  and  95 
percent  confidence  intervals  are  plotted  on  Figure  16  as  log  strength 
vs.  log  of  fiber  length.  Two  groups  of  specimens  were  tested  at 
both  the  3 and  10  in.  gage  lengths.  For  constant  diameter  fibers, 
the  volume  (and  the  area)  is  proportional  to  the  length,  therefore 
the  abscissa  of  Figure  16  is  directly  proportional  to  a log  volume 
scale.  For  a material  with  a single  distribution  of  flaws  obeying 
two-parameter  Weibull  statistics,  average  strengths  should  fall  on 


Avco,  Space  Systems  Division,  Lowell,  Mass. 


a straight  line  (with  a slope  of  -1/m)  when  plotted  in  this  manner. 
Instead,  the  data  of  Figure  16  appear  to  fall  on  two  intersecting 
straight  lines  as  expected  for  a bimodal  distribution  of  flaws. 

The  parameters  for  the  two  straight  lines  tentatively  drawn  on 
Figure  16  are  m^  = 35,  m2  = 4.0,  = 524  kpsi,  and  0^2  “ 50.3  kpsi. 

Unfortunately,  when  average  strengths  are  plotted  in  the 
manner  of  Figure  16,  much  important  information  is  discarded  concerning 
the  scatter  of  the  data  about  each  of  the  averages.  Through  the 
use  of  Equation  23,  every  data  point  can  be  visualized  simultaneously, 
allowing  a better  judgement  of  the  applicability  of  weakest  link 
statistics,  two-parameter  Weibull,  single  or  multiple  flaw  distributions, 
etc.  The  300  data  points  of  Kotchick  et  al.  were  normalized  to  a 
one  in.  gage  length  and  all  plotted  on  Figure  17  as  In  In  — 

In  strength.  The  two  data  sets  at  both  3 and  10  in  lengths  were 
combined  before  normalizing.  The  solid  line  included  in  Figure  17 
was  plotted  according  to  Equation  20  using  the  same  m's  and  a^'s 
listed  above  and  derived  from  Figure  16.  In  some  regions  of  Figure  17 
the  data  points  are  so  densely  packed  that  individual  symbols  are 
difficult  to  identify,  but  such  a representation  is  still  useful  to 
concisely  visualize  the  trends  of  the  data.  For  instance,  when 
plotted  in  this  manner,  it  appears  that  the  distribution  with  the 
lower  modulus  is  not  adequately  described  by  an  m of  4.0  as  was  estimated 
from  Figure  16.  A distribution  with  an  m of  approximately  5 would 
be  more  suitable. 

If  the  non-linear  trend  of  the  data  in  Figure  17  is  due  to  two 
active  flaw  distributions,  then  fracture  origins  of  specimens 
failing  in  one  region  of  the  curve  should  be  predominated  by  one 
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Figure  16.  Fracture  Strength  as  Measured 
in  Tension  Versus  Gage  Length 
for  Fibers  of  CVD  SiC  as  Reported 
by  Kotchick  et  al. 


Ln  STRENGTH  (ksi) 

Figure  17.  Weibxall  Plot  for  CVD  SiC  Fibers  Where  All  Probabil 
ities  Have  Been  Normalized  to  That  for  Specimens 
with  One  Inch  Gage  Lengths.  Original  Data  from 
Kotchik  et  al. 
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type  of  flaw  and  specimens  failing  in  the  other  region  should 
originate  from  a different  type  of  flaw.  Unfortunately,  this 
complimentary  data  on  fracture  origins  is  not  available  for  the  SiC 
fibers,  therefore  the  presence  of  two  active  flaw  distributions  is 
not  confirmed  and  the  non-linear  trend  of  the  data  on  Figure  17 
may  be  due  to  other  factors. 

11  12 

The  work  on  AI2OJ  by  Bansal,  Duckworth  and  Niesz  ’ involved 
a total  of  76  bend  specimens  of  Alsimag  614*  in  six  different 
variations  of  specimen  size,  load  configuration  and  environment. 
Specimens  were  prepared  in  two  sizes:  0.1  x 0.2  x 1.5  in.  and  0.5 
X 1.0  X 7.5  in.  The  smaller  specimens  were  tested  in  both  three- 
point  bending  on  a 1.25  in.  span  and  four-point  bending  with  a 1.25 
in.  outer  span  and  0.75  in.  inner  span.  The  larger  specimens  were 
all  tested  in  four -point  bending  with  a 6.25  in.  outer  span  and  a 
3.75  in.  inner  span.  The  above  specimen  configurations  were  tested 
in  two  different  environments.  Dry  nitrogen  was  used  at  high 
loading  rates  (100  MPa/s)  to  minimize  subcritical  crack  growth,  and 
distilled  water  was  used  at  slower  loading  rates  (4  MPa/s)  to 
enhance  the  rate  of  crack  growth.  The  results  of  these  six  variations, 
as  included  in  Reference  12 , are  plotted  on  Figure  18  as  log  log  "J^p" 
log  strength  along  with  the  Weibull  modulus  for  each  group  of  specimens. 
As  expected,  the  specimens  tested  in  water  had  lower  strengths  than 
those  tested  dry.  Also,  the  larger  specimens  consistently  yielded 
lower  strengths  than  small  specimens.  The  Weibull  moduli  of  the 
three  groups  tested  in  water  were  constant  at  a value  of  34,  but 
the  moduli  of  the  specimens  tested  dry  varied  from  34  for  the  small 
* 

3M  Co.,  Technical  Ceramic  Products  Division,  Chattanooga,  Tenn. 
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three-point  specimens  to  11  for  the  large  four-point  specimens. 

Fractography  as  reported  by  Bansal  et  al.  revealed  that  all 
specimens  tested  in  water  failed  from  flaws  extending  from  the 
tensile  surface.  No  voids  or  inclusions  were  associated  with  these 
flaws.  Approximately  80  percent  of  the  small  specimens  tested  in 
dry  nitrogen  also  failed  from  origins  extending  from  the  tensile 
surface,  y^proximately  10  percent  of  the  origins  were  located  at 
large,  subsurface  pores,  and  the  remaining  10  percent  could  not  be 
clearly  identified.  All  of  the  large  specimens  tested  in  dry 
nitrogen  failed  from  large,  irregular,  subsurface  pores.  Bansal  et 
al.  concluded  from  the  decrease  in  Weibull  modulus  and  the  change 
in  predominating  flaw  type  with  increasing  effective  voliime  when 
tested  in  dry  nitrogen  that  there  were  two  active  flaw  distributions. 

No  attempts  were  made  to  determine  if  the  data  were  consistent  with 
relationships  such  as  Equation  21. 

Through  the  use  of  Equation  24,  all  of  the  fracture  data  of  Figure 
18  were  normalized  to  probabilities  of  failure  for  small  three- 
point  bend  specimens  and  were  re-plotted  on  Figure  19  as  In  In  vs. 

In  fracture  strength.  Using  the  hypothesis  that  both  area  and 
volume  flaws  were  active  in  dry  nitrogen,  the  four  Weibull  parameters 
of  the  bimodal  flaw  distribution  were  needed  before  the  data  could 
be  normalized.  The  parameters  were  varied  iteratively  to  find  that 
combination  with  the  least  scatter  between  the  measured  and  predicted 
behaviors.  For  the  specimens  tested  in  dry  nitrogen,  these  parameters 
were  m^  ■ 29.1,  = 51.2  kpsi,  m2  = 8.2,  and  0^2  “ kpsi. 

The  first  two  parameters  describe  the  surface  or  area  distribution 
of  flaws,  while  the  second  two  describe  the  volume  related  flaws. 
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Figure  18.  Weibull  Plot  for  AI2O3  Showing  Six  Different  Combina- 
tions of  Specimen  Size,  Testing  Configuration  and  Test- 
ing Environment  as  Reported  by  Bansal  et  al. 
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The  solid  curve  drawn  through  the  dry  nitrogen  data  of  Figure  19 
represents  the  predicted  behavior  of  the  combined  distribution 
according  to  the  above  parameters.  This  curve  describes  the  data 
very  well  and,  when  combined  with  the  fractography  results,  provides 
a very  convincing  argument  that  the  material  contains  two  active 
flaw  distributions  and  that  each  distribution  can  be  described 
satisfactorily  by  two-parameter  Weibull  statistics. 

TTie  same  two  distributions  of  flaws  are  present  in  specimens 
tested  in  water.  If  a mechanism  of  subcritical  crack  growth  becomes 
active  due  to  the  environment,  however,  the  effects  should  be  seen 
only  on  those  flaws  on  the  surface  since  internal  flaws  are  isolated 
from  the  external  environment  in  a dense  material.  The  fracture 
origins  of  specimens  tested  in  water  were  reported  to  be  exclusively 
surface  flaws  with  no  failures  observed  at  internal  voids.  All  of 
the  fracture  data  of  specimens  tested  in  water  are  consistent  with 
a single  distribution  of  surface  flaws.  Presumably,  these  flaws 
are  of  the  same  distribution  as  the  surface  flaws  observed  in  dry 
nitrogen,  but  in  the  aqueous  environment  the  flaws  were  significantly 
extended  by  subcritical  crack  growth  prior  to  failure.  The  curve 
drawn  through  the  solid  data  points  of  Figure  19  was  calculated  by 
assuming  that  the  volume  flaws  observed  in  dry  nitrogen  were  also 
present  in  the  specimens  tested  in  water  and  that  they  were  not 
affected  by  the  aqueous  environment.  The  two  parameters  for  the 
volume  distribution  were  therefore  assumed  to  be  identical  to  those 
listed  above  for  dry  nitrogen.  The  remaining  two  parameters  for 
the  surface  distribution  were  varied  iteratively  to  find  those 


= 37.9  kpsi,  m2  = 8.2  and 


= 35.2  kpsi.  As  seen 


on  Figure  19,  the  volume  distribution  of  flaws  should  not  become 
active  within  the  range  of  probabilities  of  failure  of  the  data. 

If  more  data  had  been  collected  on  the  large  four-point  specimens 
tested  in  water  or  if  even  larger  specimens  were  tested,  then  one 
would  expect  that  some  failures  would  originate  from  the  type  of 
volume  flaws  that  were  observed  in  the  large  four-point  specimens 
tested  in  dry  nitrogen. 

In  summary,  a method  was  described  to  simultaneously  illustrate 
fracture  data  from  more  than  one  specimen  size  and  geometry  on  a 
single  two-dimensional  plot.  This  technique  proves  to  be  particularly 
useful  when  a material  contains  multiple  independent  flaw  distributions 
such  as  boron-doped  sintered  SiC.  Due  to  insufficient  data  on  sintered 
SiC  at  present,  however,  the  method  was  illustrated  quantitatively  on 
two  sets  of  literature  data  which  appear  to  contain  multiple  flaw 
distributions. 
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VII.  TIME-TEMPERATURE  DEPENDENCE  OF  FRACTl'RE  STRENGTC 

One  of  the  most  important  properties  of  a potential  high 
temperature  structural  ceramic  is  the  degree  of  strength  degradation 
which  occurs  due  to  subcritical  crack  growth  at  high  temperatures. 

All  materials  undergo  subcritical  crack  growth  to  some  degree  at 
high  temperature,  thus  resulting  in  a time-dependence  of  fracture 
strength.  The  velocity  of  subcritical  crack  growth  for  a given 
material  is  a complex  function  of  temperature,  stress  and  environ- 
ment (gas  phase  composition) . At  constant  temperature  and  environ- 
ment, the  velocity  of  crack  growth,  v,  is  often  described  empirically 

13  14 

as  a power  function  of  the  stress  intensity  factor,  K ’ 

V = A k"  (26) 

where  A and  n are  experimentally  determined  material  constants. 

Crack  velocity  can  be  measured  directly  in  tests  such  as  the 
double  cantilever  beam  and  the  double  torsion  tests,  however,  ex- 
perimental problems  can  arise  in  these  direct  measurements  at  very 
low  crack  velocities.  Indirect  methods  such  as  constant  stress 
(stress  rupture)  and  constant  loading  rate  tests  can  also  be  used 
to  evaluate  the  parameters  of  Equation  26  and  are  more  amenable  to 
testing  at  very  low  crack  velocities  (long  times  to  failure) . 

It  has  been  shown  elsewhere^^  that  Equation  26  can  be  integrated 
over  the  total  time  to  failure  at  a constant  stress  thus  yielding 
the  relationship: 

a"T  = k (27) 

where  o is  the  applied  stress,  t is  the  time  to  failure  and  k is  a 
constant  which  includes  the  critical  stress  intensity  factor,  the 
"quick-time"  fracture  strength  and  the  constant.  A,  of  Equation  26. 
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In  a similar  manner.  Equation  26  can  be  integrated  over  the  total 
time  to  failure  in  a constant  loading  rate  test  (conventional  mode 
of  bend  testing),  thus  yielding  a similar  relationship: 

o^T  = k(n+l)  (28) 

According  to  Equations  27  and  28,  stress  rupture  and  constant 
loading  rate  data  should  fall  on  straight  lines  when  plotted  as  log 
fracture  stress  versus  log  time  to  failure  as  shown  schematically 
on  Fig.  20.  The  linear  slope  in  both  cases  should  equal  - ^ and 
the  horizontal  separation  of  the  two  lines  should  be  log  (n+1) . 

Plotting  the  data  in  this  manner  provides  a simple  graphical  means 
of  predicting  the  stress  rupture  behavior  from  constant  stress-rate 
data,  but  the  relationships  in  Figure  20  are  subject  to  the  validity 
of  Equation  26. 

Charles^^  has  recently  proposed  an  expression  for  crack  velocity 
which  was  derived  in  a more  fundamental  manner  than  Equation  26 
from  a model  of  stress-directed  vacancy  diffusion  at  crack  tips. 

The  crack  velocity  predicted  from  this  model  can  be  simplified  to  the 
following  form: 

V = ^ exp  8 K (29) 

where  x is  the  instantaneous  crack  length  and  A and  8 are  material 
constants  which  can  be  related  to  thermodynamic  properties  through 
the  derivation.  The  unique  aspect  of  this  expression  is  the  crack 
length  in  the  denominator  of  the  pre-exponential.  The  model  assumes 
that  the  rate  of  crack  growth  is  controlled  by  diffusion  of  vacancies 
from  the  bulk  of  the  specimen  to  the  crack  tip  (or  conversely,  by 
diffusion  of  atomic  species  from  the  crack  tip  into  the  specimen) . 

The  crack  length  in  the  denominator  arises  because  the  average  diffusion 
distance  is  believed  to  vary  in  direct  proportion  to  the  crack  length. 


(If  Equation  29  is  correct,  it  is  interesting  to  note  that  stress 
intensity  versus  crack  velocity  curves  which  are  common  in  the 
literature  require  a third  axis  of  crack  length). 

Equation  29  has  been  integrated  for  the  case  of  constant 

stress^^.  TTie  resulting  relationship,  when  plotted  on  the  axes  of 

Figure  20,  is  concave  downward  with  an  initial  slope  (short  times 

to  failure)  that  is  a function  of  the  thermodynamic  properties  of 

the  material  and  a final  slope  (at  very  long  times  to  failure)  of  - 

1/4.  The  inverse  fourth  power  dependence  at  long  times  to  failure 

is  predicted  to  be  common  to  all  materials  — metals  and  ceramics 

alike.  Unlike  the  case  of  constant  stress.  Equation  29  cannot  be 

integrated  for  the  case  of  linearly  increasing  stress  without 

resorting  to  numerical  integration.  The  resulting  relationship, 

when  plotted  on  the  axes  of  Figure  20,  is  also  concave  downward  and 

shifted  to  the  right  of  the  constant  stress  behavior,  but  the 

amount  of  shift  is  not  simply  log  (n+1)  as  it  was  for  the  relationships 

derived  from  Equation  26.  Using  data  available  from  the  literature, 

Charles  has  shown  that  Equation  29  is  remarkably  accurate  in  describing 
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the  time- temperature  dependences  of  strength  in  both  metals  ’ 

. 17,18 

and  ceramics 

2 

Results  included  in  last  year's  final  report  indicated  that 
boron-doped  sintered  SiC  has  excellent  resistance  to  high  temperature 
subcritical  crack  growth  when  compared  with  other  potential  structural 
ceramics.  This  area  of  strength  degradation  at  high  temperatures 
is  sufficiently  inq)ortant  for  the  eventual  application  of  sintered 
SiC,  however,  that  work  has  been  expanded  during  the  previous  year. 

The  areas  to  be  discussed  in  this  report  include: 

• The  stress  rupture  behavior  of  pressurized  tubes  at 
temperatures  in  the  vicinity  of  2100®C. 
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• The  temperature  dependence  of  fracture  strength  as 
measured  in  bending  from  room  tenperature  to  IbOO^C. 

• The  dependence  of  fracture  strength  on  the  rate  of 
loading  in  three-point  bending  at  1600®C  with  times  to 
failure  ranging  from  one  second  to  five  hours. 

In  the  first  of  these  three  areas,  slip  cast  and  sintered 
tubes  of  Sic  were  tested  at  extremely  high  temperature  in  biaxial 
tension  by  locating  the  tubes  in  a graphite  resistance  furnace  and 
pressurizing  them  with  nitrogen.  The  resulting  stress  rupture  data 
along  with  an  analysis  of  the  data  have  been  written  as  an  internal 
report  which  is  included  here  as  Appendix  A.  The  data  analysis 
employs  the  stress  rupture  theory  of  Charles,  therefore  another 
internal  report  describing  the  theory  in  more  detail  is  included  as 
^pendix  B. 

The  second  and  third  areas  of  time-temperature  dependent  strength 
listed  above  were  performed  on  a common  set  of  approximately  90  bend 
specimens  which  were  prepared  in  a manner  similar  to  that  described 
in  Section  II.  The  bars  were  sintered  in  vacuum  at  a temperature  of 
approximately  2080®C  for  30  minutes,  resulting  in  densities  from  95 
to  96  percent  of  theoretical  and  microstructures  conposed  primarily 
of  3-5pm  grains  of  g-SlC.  The  bars  were  carefully  ground  from  the 
as-sintered  cross-sectional  dimensions  of  0.167.x  0.167  in.  to 
0.150  X 0.150  in.  with  a slight  chamfer  on  the  edges.  The  bars 
were  not  refired  before  testing.  All  testing  (including  room 
temperature)  was  carried  out  in  three-point  bending  on  a one  and 
one-half  inch  span  using  a carefully  aligned  fixed-pin  test  fixture 
fabricated  from  sintered  SiC.  An  internally  wound  Pt-40%  Rh  tube 
furnace  allowed  testing  in  air  at  temperatures  to  1600°C.  The 
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furnace  and  fixtures  were  mounted  in  an  Instron  testing  machine 
capable  of  crosshead  speeds  from  0.00002  to  20  in. /min. 

The  ten^jerature  dependence  of  strength  was  studied  by  testing 
ten  specimens  at  each  of  four  temperatures  from  room  temperature  to 
1600“C  at  a crosshead  speed  of  0.002  in. /min.  and  a corresponding 
loading  rate  of  approximately  17  kspi/min  (115  MPa/min) . The 
resulting  strengths  are  plotted  on  Figure  21  as  a function  of 
temperature.  Again,  the  data  points  represent  the  average  strengths 
while  the  error  bars  include  the  95  percent  confidence  limits  on 
the  averages.  Hie  strength  appears  to  be  independent  of  temperature 
dt  a stress  of  approximately  80  kpsi  (550  MPa)  from  room  temperature 
to  about  1200®C  and  drops  to  approximately  70  kpsi  (480  MPa)  by 
1600®C.  Examination  of  the  fracture  surfaces  indicated  that  the 
fracture  origins  were  located  at  pro cessing -related  flaws  in  the 
microstructure  rather  than  machining  damage. 

The  loading  rate  dependence  of  fracture  strength  in  sintered 
SiC  was  studied  in  air  at  1600°C  in  order  to  estimate  the  parameters 
of  crack  velocity  relationships  such  as  Equations  26  and  29. 

Constant  loading  rate  tests  were  chosen  rather  than  constant  stress 
tests  for  the  majority  of  the  present  work  because  of  experimental 
constraints  of  the  testing  apparatus.  The  testing  procedure  included: 
heating  the  furnace  to  1600®C,  raising  the  furnace  and  upper  fixture 
to  expose  the  lower  testing  fixture,  positioning  the  specimen  on 
the  lower  fixture,  lowering  the  furnace  and  upper  fixture  into  the 
testing  position,  equilibrating  the  temperature  for  15  minutes,  and 
loading  at  a constant  crosshead  speed  until  failure.  Five  crosshead 
speeds  were  chosen  from  0.5  to  0.00005  in. /min  resulting  in  times 
to  failure  of  one  second  to  five  hours.  At  least  five  specimens 
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were  tested  at  each  loading  rate.  The  results  are  included  as  the 
open  circles  on  Figure  22  where  log  fracture  strength  is  plotted  as 
a function  of  log  time  to  failure  (same  axes  as  Figure  20) . The 
fastest  loading  rates  resulted  in  an  average  strength  of  82.5  kpsi 
(570  MPa)  which  is  effectively  the  same  as  the  room  ten5)erature 
strength  included  on  Figure  21.  As  expected,  slower  loading  rates 
and  longer  times  to  failure  resulted  in  steadily  decreasing  fracture 
strengths.  The  data  points  fall  on  a smooth  curve  which  is  concave 
downward  as  predicted  by  Equation  29  instead  of  linear  as  predicted 
by  Equation  26.  Unfortunately,  the  rate  of  decrease  in  the  slope 
of  the  curve  as  drawn  on  Figure  22  is  greater  than  can  be  explained 
conveniently  by  Equation  29.  Also,  it  was  observed  that  the  specimens 
tested  at  the  slowest  crosshead  speed  of  0.00005  in. /min.  were 
covered  with  a thick  layer  of  spalling  oxide  which  contained  metallic 
inclusions.  The  metallic  inclusions  were  identified  by  X-ray 
diffraction  as  a Pt-Rh  alloy  which  had  probably  evaporated  from  the 
internal  windings  of  the  furnace  and  condensed  on  the  cooler 
specimen  (similar  metallic  particles  also  covered  those  regions  of  the 
test  fixtures  located  in  the  hot  zone) . 

Testing  continued  with  six  constant  stress  tests  at  the  average 
failure  stress  of  the  0.00005  in. /min.  specimens  (45.8  kpsi  - 315  MPa). 
These  results  are  included  on  Figure  22  as  the  single  square  data  point 
with  the  horizontal  error  bar.  The  average  time  to  failure  was  sur- 
prisingly long  at  approximately  75  percent  of  the  constant  loading 
rate  tests.  Relationships  such  as  Equations  26  and  29  predict  that 
the  time  to  failure  in  constant  stress  should  fall  between  2 and  10 
percent  of  the  time  to  failure  in  constant  loading  rate.  Again,  the 
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specimen  surfaces  were  severely  oxidized. 

The  oxidation  observed  in  four  to  five  hours  at  ISOO^C  in  the 

Pt-40%  Rh  furnace  is  considerably  more  severe  than  similar  times 
and  temperatures  in  other  furnaces.  This  aggrevated  oxidation  and 
spalling  may  be  the  result  of  contamination  of  the  protective  Si02 
layer  with  a platinum  silicide  of  some  type  or  with  AI2O2  from  the 
fibrous  insulation  used  on  the  ends  of  the  furnace.  In  any  case, 
the  oxidation  degraded  the  specimen  surfaces  sufficiently  to  suggest 
that  oxidation  contributed  to  the  strength  degradation  of  all 
specimens  tested  at  long  times  to  failure.  To  test  this  hypothesis, 
a set  of  five  specimens  were  tested  at  a crosshead  speed  of  0.020  in. /min. 
in  a similar  manner  to  the  earlier  group  tested  at  the  same  loading 
rate.  In  this  case,  however,  the  specimens  were  allowed  to  soak  in 
the  furnace  for  five  hours  before  loading.  TTie  results  are  plotted  on 
Figure  22  as  the  single  triangular  data  point.  TTie  average  strength  dropped 
from  77.3  kpsi  (530  MPa)  to  60.3  kpsi  (415  MPa)  due  to  the  additional  five 
hours  of  exposure  in  the  furnace. 

The  data  discussed  above  confirm  that  two  mechanisms  are  active  in 
decreasing  the  high  temperature  fracture  strength  with  increasing  testing 
time.  The  first  is  a mechanism  of  subcritical  crack  growth  and  the 
second  is  a surface  degradation  due  to  oxidation.  It  is  difficult 
to  separate  the  two  mechanisms  with  the  data  available,  but  two  related 
attempts  can  be  made  to  salvage  some  information  concerning  subcritical 
crack  growth.  First,  the  dashed  straight  line  of  Figure  22  connects 
two  constant  loading  rate  data  points  which  had  similar  histories  of 
approximately  five  hours  in  the  furnace  at  1600®C  prior  to  failure. 

From  Equation  26,  the  slope  of  this  line  yields  an  estimate  of  the 
exponent,  n,  of  approximately  25.  The  second  method  is  somewhat  less 
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Straightforward.  An  additional  five  hours  of  furnace  exposure 

prior  to  failure  at  a crosshead  speed  of  0.020  in. /min.  decreased 

the  average  strength  by  approximately  22  percent.  Conversely,  the 

hypothetical  case  of  testing  at  0.00005  in. /min.  without  the  five 

hours  of  furnace  exposure  during  the  test  should  increase  the  strength 

by  a comparable  amount.  The  "X"  included  on  Figure  22  represents  this 

hypothetical  data  point  which  was  tested  at  0.00005  in. /min.  and  had 

a time  to  failure  of  approximately  five  hours  but  was  located  in  the 

testing  furnace  for  only  20  to  30  minute  . This  is  the  approximate 

time  of  furnace  exposure  for  the  remaining  four  similar  data  points 

on  Figure  22  (with  the  exception  of  the  0.0005  in. /min.  specimens 

which  were  in  the  furnace  for  approximately  45  minutes) . A least- 

squares  straight  line  through  these  five  points  (not  drawn  on  Figure 

22)  yields  an  exponent,  n,  of  approximately  28.  Both  of  these  estimates 

2 

are  reasonably  consistent  with  earlier  measurements  of  33. 

Table  1 includes  all  of  the  average  strengths,  standard 
deviations,  times  to  failure,  etc.  for  the  data  reported  on  Figures 
21  and  22. 

In  summary,  the  time- temperature  dependence  of  fracture  strength 
in  boron-doped  sintered  SiC  was  measured  in  three-point  bending.  The 
strength  was  found  to  be  independent  of  temperature  from  room  tem- 
perature to  approximately  1200®C  and  dropped  to  approximately  85  per- 
cent of  the  room  temperature  strength  by  1600®C.  The  time  dependence 
of  strength  at  1600®C  was  studied  with  constant  loading  rate  tests. 

In  addition  to  subcritical  crack  growth,  a second  mechanism  of  high 
tenq>erature  strength  degradation  was  active  and  appeared  to  be  related 
to  oxidation.  This  second  mechanism  is  believed  to  be  characteristic 
of  the  particular  furnace  used  during  this  testing.  Similar  conditions 
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TABLE  1 


RESULTS  OF  HIGH  TEMPERATURE  BEND 
TESTING  USING  CONSTANT  LOADING  RATES 


Temp. 

(’C) 

Crosshead 
Speed 
(in. /min.) 

Avg.  Time  to 
Failure 
(sec.) 

Avg.  Strength 
(kpsi) 

Std.  Deviation 
(kpsi) 

Number  of 
Specimens 

RT 

0.010 

50 

80.4 

7.3 

10 

1200 

0.002 

250 

77.4 

10.2 

11 

1400 

0.002 

270 

76.7 

7.2 

10 

1600 

0.002 

360 

70.3 

3.9 

10 

1600 

0.500 

1.2 

82.5 

17.4 

10 

1600 

0.020 

33.8 

77.3 

7.6 

5 

1600 

0.0005 

1510 

62.5 

6.5 

5 

1600 

0.00005 

22,000 

45.8 

2.9 

10 

1600* 

— 

16,800 

45.8 

— 

6 

lOOO"^ 

0.020 

30.2 

60.3 

9.3 

5 

I 


Constant  stress 

^Loading  began  after  five  hours  of  soak  in  the  testing  apparatus 
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in  other  furnaces  lead  to  markedly  less  surface  degradation  due  to 
oxidation  and  spalling. 
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VIII.  CREEP  MEASUREMENTS 


I 


Creep  of  boron-doped  silicon  carbide  was  studied  previously  on 

hot-pressed  specimens  where  strain  rates  were  measured  in  three- 

19 

point  bending  near  1600®C  in  air  at  a stress  of  25  kpsi.  Creep 

-9  -8 

rates  were  measured  between  1.5  x 10  and  1 x 10  /sec.  with 
stress  exponents  between  1.2  and  1.6.  Subsequent  measurements  of 
fine-grained  sintered  SiC  of  about  95  percent  theoretical  density 
yielded  no  steady-state  creep  rate  even  after  150  hours  at  1550®C. 
Because  of  this  uncertainty  the  results  were  not  reported. 

The  closed-end-tube  specimen  geometry  adopted  for  stress - 
n:5)ture  tests  described  in  ^pendix  A offers  the  possibility  of 
creep  measurements  under  those  conditions  that  cause  measurable 
strain  (change  in  diameter  of  the  tube)  before  rupture.  This  was 
not  the  case  with  the  specimens  used  for  the  rupture  tests  described 
in  i^pendix  A where  the  tubes  were  sintered  at  2150‘’C,  were  coarse 
grained,  and  showed  less  than  0.010  in.  of  radial  deformation  after 
exposure  to  an  internal  pressure  of  1350  psi  (9.3  MPa)  at  2100°C 
for  45  minutes.  Therefore  another  set  of  tubes  was  slip  cast  and 
sintered  at  2080“C  for  40  minutes  in  vacuum.  This  temperature  was 
selected  as  the  lowest  at  which  helium  tight  tubes  were  reliably 
obtained.  Although  a high  degree  of  3 to  a transformation  occurred, 
a fairly  uniform  microstructure  with  no  "feathers"  resulted  as 
shown  in  Figure  23. 

For  creep  measurements  the  tubes  were  bonded  with  epoxy  into 
the  pressure  head,  proof-tested  with  a gas  pressure  of  2200  psi 
(15.2  MPa)  at  room  temperature  and  inserted  into  the  furnace.  The  tubes 
were  positioned  in  the  furnace  so  that  the  hottest  region  was 
about  0.8  in.  from  the  closed  end.  This  established  a 50®C  difference 
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between  the  temperature  of  the  hot  zone  and  the  temperature  measured 
at  the  closed  end  and  decreased  the  possibility  of  any  end  effects 
in  the  measurement  of  strain.  Temperature  was  measured  as  described 
in  Appendix  A using  a Leeds  amf  Northriq)  optical  pyrometer  which 
was  calibrated  at  the  melting  point  of  sapphire  in  argon  in  the 
same  furnace.  The  surface  carbon  (from  sintering  in  vacuum)  was 
oxidized  from  each  specimen  with  a torch  and  the  surface  was  lightly 
sanded  with  600  grit  abrasive  paper.  On  each  run  the  tube  was 
inserted  into  the  furnace,  brc'ight  to  tenqjerature  in  about  10 
minutes,  equilibrated  for  another  10  minutes,  and  pressurized  with 
nitrogen  to  1200  psi  (8.3  MPa)  for  the  desired  time  interval  (15 
to  30  minutes).  After  depressurizing  the  tube,  the  furnace  was 
cooled,  the  tube  was  removed  and  its  diameter  at  the  hot  zone  was 
measured  with  a dial  gauge  to  within  0.0001  in.  The  tube  was 
slowly  rotated  through  360  degrees  in  a V-block  and  readings  were 
taken  at  the  largest  and  smallest  diameter.  The  resulting  four 
readings  were  averaged. 

After  collecting  sufficient  data,  a 2.5  in.  portion  of  the 
exposed  end  of  the  tube  was  removed  and  cut  longitudinally  into  two 
halves.  One  half  was  polished  and  etched  for  metallography  while 
the  other  was  sliced  into  1/4  inch  long  segments  which  were  evaluated 
by  liquid  displacement  for  a density  profile. 

Three  specimens  have  been  tested  to  date,  all  at  1200  psi  (8.3  MPa) 
internal  nitrogen  pressure  with  an  external  atmosphere  of  nitrogen 
at  ambient  pressure.  The  results,  in  terms  of  change  in  diameter  and 
strain  are  included  in  Table  2.  Table  3 includes  initial  dimensions, 
wall  thicknesses,  and  wall  stresses  of  the  specimens.  The  results  are 
also  plotted  on  Figure  24  as  strain  versus  time. 

Specimen  No.  1 was  tested  at  2010®C  (curve  1,  Figure  24). 
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TABLE  2 

I 

SUMMARY 

OF  CREEP  EXPERIMENTS 

t 

Exp.  # 

Temp. 

("CD 

Time  at  Ten^.  and 
Press.  (Min.) 

Cumulative  Time 
(Min.) 

Change  in  Oia. 
(in.  X 10 3) 

Strain 

- 

TUBE  1 

1 

2010 

30 

30 

1.5 

.0024 

2 

2010 

30 

60 

2.2 

.0043 

1 

i 

3 

2010 

30 

90 

2.6 

.0050 

( 

4 

2010 

30 

120 

3.15 

.0061 

5 

2010 

30 

150 

3.65 

.0070 

6 

2010 

30 

180 

4.25 

.0082 

7 

2010 

30 

210 

4.85 

.0094 

8 

2010 

30 

240 

5.55 

.0107 

9 

2170 

30 

— 

— 

10 

2070 

30 

270 

5.55 

.0107 

1 

11 

2120 

30 

300 

5.75 

.0111 

12 

2120 

30 

330 

6.00 

.0115 

13 

2120 

30 

360 

6.15 

.0117 

- : 

TUBE  2 

I 

14 

2070 

15 

15 

3.65 

.0072 

15 

2070 

15 

30 

4.35 

.0087 

16 

2070 

15 

45 

5.10 

.0103 

17 

2070 

15 

60 

6.35 

.0126 

18 

2070 

15 

75 

7.75 

.0148 

19 

2070 

15 

90 

8.54 

.0169 

20 

2070+ 

15 

105 

10.43 

.0205 

21 

2C70 

15 

120 

11.50 

.0226 

22 

2070 

15 

135 

12.40 

.0243 

23 

2070 

20 

155 

13.85 

.0272 

24 

2070 

15 

170 

15.80 

.0306 

25 

2070 

15 

185 

18.40 

.0359 

TUBE  3 

26 

2070 

45 

45 

5.82 

.0113 

27 

2120 

15 

60 

9.80 

.0190 

28 

2120 

15 

75 

14.95 

.0295 

' 

29 

2120 

15 

90 

21.65 

.0419 
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CREEP  STR/ 


Figure  24,  Creep  Strain  Versus  Time  for  Pressurized  Tubes 
at  Three  Different  Temperatures 


After  an  initial  60  minute  period  of  primary  creep  (decreasing 

creep  rate),  the  data  points  seem  to  fall  on  a straight  line  which 

-4  -8 

yields  a steady-state  creep  rate  of  0.36  x 10  /min.  (1  x 10  /sec.). 
After  4 hours  of  exposure  at  2010°C  the  tube  was  heated  in  argon 
(not  pressurized)  to  2170‘’C  for  30  minutes  to  bring  about  substantial 
grain  growth.  During  this  exposure  some  evaporation  of  the  outer 
surface  took  place  so  that  the  diameter  actually  decreased  by  0.005 
in.  This  change  in  diameter  has  been  taken  into  account  in  data 
points  10-13  of  Table  2.  The  tube  was  then  pressurized  at  2070®C 
for  30  minutes  with  no  detectable  deformation.  The  temperature  was 
increased  to  2120®C  where  some  creep  could  be  detected.  It  was 
realized,  however,  that  these  very  slow  creep  rates  would  be  affected 
by  simultaneous  evaporation,  thus  the  experiment  was  discontinued. 

The  microstructure  of  this  specimen  sectioned  in  the  hot  zone  is 
shown  in  Figure  25.  Comparison  with  the  initial  microstructure  of 
Figure  23  suggests  a strong  grain  size  dependence  of  creep  rate  in 
boron-doped  SiC. 

Specimen  2 was  tested  at  2070*0  in  15  minute  intervals  for  3 
hours.  Curve  2 on  Figure  24  shows  the  strain  versus  time  behavior 
and  indicates  steady- state  creep  after  only  15  minutes  of  testing. 

During  the  seventh  15  minute  pressurization  the  temperature  accidentally 
increased  about  25 ®C  for  a short  but  unknown  period  of  time  which 
resulted  in  accelerated  creep  and  a discontinuity  in  the  linear 
behavior  as  plotted  on  Figure  24.  At  about  3 percent  strain  (160 
minutes  total  exposure)  the  strain  rate  began  increasing,  thus 
indicating  the  onset  of  tertiary  creep.  From  the  linear  segment  of 
curve  2 one  obtains  a steady-state  creep  rate  at  2070®C  of  1.38  x 
10"Vrain.  (3.83  x 10"®/sec) . After  termination  of  the  experiment. 
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the  tube  was  sectioned,  polished  and  etched  for  optical  microscopy 
and  evaluation  of  density.  Figure  26  shows  an  etched  section  of 
the  tube  exposed  to  the  highest  temperature.  No  obvious  grain 
growth  occurred  (compare  with  Figure  23)  as  a result  of  the  testing. 

The  longitudinal  density  profile  is  shown  on  Figure  27  as  porosity 
versus  distance  from  the  end  of  the  specimen  showing  a one  percent 
increase  in  porosity  due  to  the  deformation  of  the  hot  zone  section. 

Specimen  3 was  first  exposed  for  45  minutes  at  2070°C.  Only 
the  final  strain  was  recorded  over  this  period  but  was  consistent 
with  data  from  Specimen  2.  The  specimen  was  then  exposed  for  three 
15  minute  periods  at  2120°C.  Creep  was  rapid  and  resulted  in  4.2 
percent  of  strain  after  45  minutes  at  which  point  the  experiment 
was  terminated.  Determination  of  a steady-state  creep  rate  from 
these  measurements  would  be  uncertain  and  has  not  been  attempted. 

The  data  are  plotted  on  curve  3 of  Figure  24.  Grain  growth  did  not 
occur  but  an  increase  in  pore  size  can  be  seen  on  as-polished 
sections  in  Figures  28a  and  b (before  and  after  testing,  respectively) 
and  also  on  the  etched  section  of  Figure  29.  Evaluation  of  density 
as  a function  of  position  revealed  an  increase  in  porosity  in  the 
hot  zone  of  about  two  percent  as  shown  on  Figure  27. 

The  measurement  of  creep  on  closed-end  tube  specimens  at  very 
high  temperature  has  the  advantage  that  it  eliminates  the  uncertainties 
related  to  the  rigidity  of  the  loading  train  and  specimen  grips 
usually  encountered  in  other  techniques.  In  addition  the  resulting 
tensile  stress  is  effectively  constant  therefore  allowing  estimates 
of  stress  and  strain  without  the  knowledge  of  the  stress  exponent. 
Assuming  a uniform  maximum  tensile  stress  of  = Pd/2t  in  the  tube 
wall  due  to  internal  pressure,  P,  and  diffusional  creep  (stress 
exponent  = 1),  one  can  calculate  the  activation  energy,  AH,  for 
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0 05  1.0  1.5  2.0 

DISTANCE  FROM  END  OF  SPECIMEN  (in.) 

Figure  27.  Density  After  Testing  Versus  Distauice  from  the  End  of  the  Tube 
for  Two  Creep  Specimens 
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(b) 

Figure  28.  As -Polished  Sections  at  the  Hot  Zone 


of  Specimen  No.  3 Illustrating  the  In- 
crease in  Porosity  from  Before  (a) 
to  After  Testing  (b).  21  OX 
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the  rate  controlling  process  as  follows: 


-)  (30) 


Substituting  the  stresses,  steady-state  strain  rates  and  test 

temperatures  for  Specimens  1 and  2 results  in  a calculated  AH  of 

246  Kcal/mole.  This  may  be  compared  to  activation  energies  of  165 

to  210  Kcal/mole  reported  previously  from  creep  data  of  hot-pressed 

19 

SiC  in  bending  near  1600®C.  The  highest  activation  energies  of 
this  earlier  study  were  measured  on  a chemical  composition  very 
similar  to  the  sintered  SiC  used  in  the  present  pressurized  tube 
experiments . 

It  is  customary  to  calculate  apparent  diffusion  coefficients 

from  creep  data  and  compare  the  results  with  diffusion  coefficients 

obtained  by  other  techniques,  particularly  by  tracer  diffusion. 

The  calculation  of  from  creep  data,  however,  requires  information 

20 

on  grain  size.  In  Herring's  treatment,  creep  in  a polycrystalline 

solid  is  modelled  by  the  flow  of  lattice  defects  which  are  motivated 

by  chemical  potential  gradients  in  the  grains.  The  grain  size 

enters  as  a squared  relationship  for  the  model  of  a uniform,  equiaxed 

grain  structure.  Raj  and  Ashby  have  discussed,^ creep  in  model 

21 

structures  conposed  of  elongated  grains.  Their  solution  for  the 
strain  rate  of  high  aspect  ratio  grain  structures  (Equation  26, 
Reference  21)  reduces  to  an  expression  very  similar  to  that  of 
Herring  if  diffusion  is  limited  to  volume  diffusion  and  grain  size 
is  replaced  by  particle  length. 

The  grain  structures  of  the  present  SiC  creep  specimens  are 
coiqmsed  of  thin  interpenetrating  randomly  oriented  plates  with  a 
wide  distribution  of  sizes.  If  we  model  such  a structure  as  an 


assembly  of  discs  of  constant  diameter  to  thickness  ratio  we  may 

22  — 

obtain,  according  to  Pullman,  the  average  diameter,  d,  of  the 

discs  from  measurements  of  the  mean  linear  intercept  and  the  area 

count  in  a section.  This  quantity  was  selected  as  the  characteristic 

dimension  for  estimation  of  D . For  thin  plates  Pullman  derived: 

app  ^ 


d = 


(31) 


Where  is  the  number  of  particles  per  unit  length  of  a section 

and  Ng  is  the  number  of  particles  per  unit  area.  This  yielded  d = 

14.6vun  for  the  sintered  tubes  of  SiC.  D can  then  be  calculated 

app 

as  follows:  ^ 


app 


0.075  e k T d ^ 


(32) 


where  0,  the  molar  volume  of  defects,  is  6.23  cm  /mole.  The  resulting 

values  of  D. 

2, 


app 


are  3.73  x lO'^^  cm^/sec.  at  2343°K  and  9.02  x lO"^^ 


cm  /sec  at  2283°K.  This  data  is  plotted  on  Figure  30  along  with 
previous  values  of  D calculated  from  creep  measurements  for  hot- 
pressed  boron-doped  material,  R-81  (Ref.  19,  p.  26).  A linear 
regression  analysis  of  the  resulting  five  data  points  yielded: 


%p  " ^ [-(225  ^ 29)/RT] 


(33) 


where  the  uncertainty  in  activation  energy  represents  the  95  percent 
confidence  limits  on  the  slope  of  Figure  30  as  calculated  using 
Student's  t-distribution. 

In  addition  to  the  creep  data,  carbon  diffusion  data  from 

single  crystals  of  n and  p-type  a-SiC  according  to  Ghoshtagore  and 
23 

Coble  have  been  included  on  Figure  30.  It  is  pertinent  to  emphasize 
that  no  a priori  agreement  between  the  two  curves  is  necessary, 
however,  since  different  rate  controlling  mechanisms  may  be  active 
in  creep  and  tracer  diffusion. 
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PRESENT  WORK  SINTERED  SiC 
PRESSURIZED  TUBE 

HOT-PRESSED  SiC  5PT.  BENDING 

Dc  IN  n-TYPE  SINGLE  CRYSTAL 

Dc  IN  p-TYPE  SINGLE  CRYSTAL 

FARNSWORTH  COBLE  CREEP 
OF  Al  DOPED  HOT-PRESSED  SiC 


l/TxlO**  (•K) 

Figure  30.  Diffusion  Coefficients  Versus  Reciprocal  Temperature  Including 
Both  Data  Calculated  from  Creep  Measurements  and  Direct  Mea 
surements  of  Diffusion  from  the  Literature 
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Inspection  of  the  pore  structure  in  the  sectioned  SiC  specimens 

before  and  after  creep  (Figure  28a  and  b)  does  not  reveal  grain 

boundary  separation  as  was  clearly  observed,  for  instance,  in  creep 

of  pure  polycrystalline  alumina  composed  of  equiaxed  grains  reported 
24 

by  Folweiler.  Such  grain  boundary  separation  is  in  general 
attributed  to  a substantial  contribution  of  non-accommodated  grain 
boundary  sliding.  In  the  present  case  the  change  in  pore  morphology 
accompanying  the  creep  of  SiC  appears  to  be  resulting  from  uniform 
pore  growth.  Whether  this  difference  in  morphology  of  the  growing 
pores  results  from  the  difference  in  grain  shape  or  from  a difference 
in  creep  mechanisms  cannot  be  resolved  from  available  information. 

It  is  believed,  however,  that  this  difference  is  substantial  and 
may  be  particularly  important  with  respect  to  the  stress  rupture 
behavior. 


. MJL. 
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ORIGIN  OF  STRENGTH  LIMITING  FLAWS 

The  optimum  processing  procedure  for  SiC  powder  presently 
includes  screening  through  an  8 um  screen  as  a dispersion  of  powder 
in  benzene,  spraying  into  liquid  nitrogen  to  freeze  the  benzene 
without  segregation  of  the  binder,  and  freeze  drying  to  remove  the 
benzene.  Great  care  is  taken  in  all  steps  after  the  screening  to 
avoid  contamination  of  the  powder  with  dust,  lint,  etc.,  but 
despite  all  efforts,  a large  proportion  of  failures  in  material 
sintered  from  such  powder  originate  at  large  irregular  voids  which 
are  commonly  20-200  pm  in  length.  Figure  31  is  a scanning  electron 
micrograph  of  a typical  fracture  origin  in  material  which  was  sintered 
from  fully  processed  powder. 

The  screening  operation  is  believed  to  remove  all  foreign  material 
and  aggregates  of  SiC  larger  than  approximately  10  pm,  therefore  only 
two  plausible  sources  of  the  flaws  remain. 

• The  powder  is  contaminated  after  the  screening  with  large 
particles  of  foreign  material  which  evaporate,  decon^ose 
or  dissolve  during  the  sintering,  thus  leaving  large  un- 


sinterable  voids  behind. 


• Agglomerates  of  SiC  (probably  very  weak  ones)  are  formed 
at  some  stage  after  the  screening  and  occasionally  cause 
massive  "bridging"  during  die-pressing.  The  large  en^ty 
cavities  or  low  density  regions  that  are  created  have  no 
chance  of  being  closed  by  sintering  and  thus  result  in 
large  voids  and  flaws. 


During  the  previous  year,  progressively  increasing  care  has  been 
taken  to  avoid  contamination  of  the  powder  with  any  foreign  material 
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after  the  screening  step.  These  efforts  do  not  appear  to  have 
decreased  the  frequency  of  occurrence  of  failure  from  large  voids, 
thus  suggesting  that  they  are  not  caused  by  contamination. 

The  design  of  an  experiment  to  unambiguously  identify  the  source 
of  these  flaws  has  proven  to  be  very  difficult.  One  experiment  which 
yielded  some  information  involved  the  hot-pressing  of  fully  processed 
powder.  It  was  reasoned  that  the  large  voids  tending  to  form  during 
sintering  would  be  closed  and  fully  densified  during  hot-pressing.  If 
these  potential  voids  were  caused  by  contamination,  the  most  likely 
candidates  would  be  carbonaceous  material  such  as  lint,  dust,  etc. 
which  would  decompose  at  relatively  low  temperatures  during  sintering 
and,  depending  on  composition,  would  leave  remnants  of  carbon  on  the 
inner  surfaces  of  the  voids.  During  hot-pressing  the  foreign  material 
would  still  decompose  and  yield  carbon,  but  the  resulting  voids  would 
tend  to  close  due  to  the  additional  driving  force  for  densification 
provided  by  the  externally  applied  pressure.  The  remnant  carbon, 
however,  would  be  trapped  as  a lamellar  layer  with  comparable  dimensions 
(in  two  directions)  to  the  original  void,  and  would  behave  as  a flaw 
of  comparable  severity  to  the  original  void.  Therefore,  if  the  flaws 
observed  in  sintered  material  are  due  to  carbonaceous  foreign  material, 
then  specimens  from  hot-pressed  bodies  should  fail  at  flaws  of  com- 
parable dimensions  and  at  comparable  strengths  to  the  sintered  material. 

On  the  other  hand,  if  the  voids  in  sintered  material  are  due  to 
"bridging"  of  weak  agglomerates  during  cold  pressing,  then  hot-pressing 
of  the  same  powder  should  close  the  voids  and  densify  the  structure 
with  no  evidence  of  the  preexisting  flaws.  Failure  should  be  caused 
by  a different  type  of  flaw  which  is  less  severe  than  the  voids  in  the 
sintered  material  and  the  fracture  strength  should  accordingly  be  higher. 
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I r A charge  of  85  grams  of  fully  processed  SiC  powder  was  hot- 

!1  pressed  in  a two  in.  diameter  graphite  die  at  18S0°C  for  60  minutes 

I at  10  kpsi  (68.9  MPa)  resulting  in  a final  density  of  greater  than 

I 

r 

I 99.5  percent  of  theoretical.  The  billet  was  cut  and  ground  into 

12  bend  specimens  with  dimensions  of  0.150  x 0.150  x 1.75  in.  with 

i 

slightly  chamfered  edges.  It  was  preferred,  of  course,  that  the 
eventual  strength  limiting  flaws  would  not  be  caused  by  the  machining, 
therefore  great  care  was  taken  to  minimize  machining  damage  and  the 
specimens  were  refired  to  reduce  the  severity  of  any  flaws  which  may 
have  been  induced.  All  12  specimens  were  refired  at  1940®C  for  40 
minutes  at  a furnace  pressure  of  0.35  torr.  The  choice  of  these  re- 
firing conditions  was  based  both  on  the  data  presented  in  Section  III 
and  on  micros tructural  data  presented  below. 

The  specimens  were  tested  at  room  temperature  in  three-point 
bending  on  a one  and  one-half  inch  span  resulting  in  an  average  strength 
of  117  kpsi  (805  MPa),  a standard  deviation  of  13.3  kpsi  (91.7  MPa),  and 
a Neibull  modulus  of  8.8.  These  can  be  con^ared  to  sintered,  machined 
and  refired  specimens  idiich  generally  yield  an  average  strength  of  80-100 
kpsi,  a standard  deviation  of  9-12  kpsi  and  a Weibull  modulus  of  6-10. 

The  fracture  origins  of  all  12  specimens  were  located  and  studied 
by  scanning  electron  microscopy.  As  hoped,  failure  was  not  predominated 
by  machining  flaws  since  only  4 of  12  specimens  failed  from  origins  on 
or  near  the  surface.  Each  of  the  remaining  specimens  appeared  to  contain 
an  inclusion  at  the  fracture  origin.  Examples  of  these  are  included  as 
scanning  electron  micrographs  in  Figures  32a  and  b.  Energy  dispersive 
X-ray  analysis  of  several  of  these  regions  revealed  high  concentrations 
of  titanium.  It  is  interesting  to  note  that  the  inclusions  at  all  8 
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(b) 

Figure  32.  Scanning  Electron  Micrographs 
of  Two  Typical  Volume  Flaws 
in  Hot-Pressed  SiC.  5000X 
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subsurface  fracture  origins  measured  approximately  S to  15  pm.  This 
is  roughly  an  order  of  magnitude  smaller  than  the  voids  commonly 
observed  as  fracture  origins  in  sintered  material  (Figure  31)  and 
is  probably  small  enough  to  pass  through  the  "nominal"  8 pm  screen 
used  in  powder  processing. 

Although  the  results  of  this  hot-pressing  experiment  are  not 
conclusive,  the  following  observations  (combined  with  the  reasoning 
discussed  above)  suggest  that  the  source  of  voids  in  sintered  material 
is  not  large  particles  of  carbonaceous  material: 

• The  average  strength  of  hot-pressed  material  was  20 
to  30  percent  higher  than  that  for  sintered  material. 

• The  fracture  origins  in  hot-pressed  material  were 
approximately  an  order  of  magnitude  smaller  than 
those  in  sintered  material. 

• Several  of  the  fracture  origins  of  hot-pressed  material 
have  been  identified  as  titanium-containing  inclusions. 

It  is  concluded  that  the  most  probable  cause  of  the  large  voids 
during  the  sintering  of  SiC  is  massive  bridging  of  agglomerates  during 
the  powder  compaction  step  of  specimen  fabrication. 

As  mentioned  above,  the  12  hot-pressed  specimens  were  refired  at 
1940®C  for  40  minutes.  The  choice  of  a refiring  temperature  for  these 
specimens  required  the  consideration  of  a factor  not  normally  required 
in  the  refiring  of  sintered  and  machined  specimens.  Since  the  hot-pressed 
material  was  consolidated  at  a maximum  temperature  of  1850°C,  the  grain 
size  was  much  finer  (<  1 pm  average)  than  that  for  material  sintered 
at  2100®C  (commonly  3-5  pm  average) . No  information  was  available  on 
the  microstructural  changes  caused  by  annealing  of  such  fine  grain  sizes 
in  boron-doped,  hot-pressed  6-SiC,  therefore  some  preliminary  experiments 
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were  carried  out  on  scraps  of  the  hot -pressed  billet. 


i 1 
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Small  samples  were  annealed  for  30  minutes  at  each  of  the 
following  temperatures:  1735,  1825,  1920,  1965,  2015,  2055,  and 
2105 °C.  The  annealed  samples  and  one  as -hot -pressed  sanple  were 
polished  and  etched  in  both  molten  KCXl  * 10%  KNO^  for  60  seconds  at 
approximately  450‘’C  and  boiling  Murakami's  reagent  for  five  minutes. 
The  samples  annealed  at  1735**  and  1825'’C  revealed  no  detectable 
change  in  microstructure  as  might  be  expected  with  a hot-pressing 
tenperature  of  1850*C.  Micrographs  of  the  as -hot -pressed  sasple 
and  the  remaining  annealed  samples  are  included  in  Figures  33a-f  in 
order  of  increasing  annealing  temperature.  Coarsening  of  the 
microstructure  is  easily  resolved  at  1965'’C  (Figure  33c)  and  proceeds 
as  a type  of  exaggerated  grain  growth  where  a few  grains  grow 
freely  at  the  expense  of  the  remaining  fine-grained  8-SiC  matrix. 

The  growing  grains  are  relatively  equiaxed  and  are  primarily  8-SiC 
as  opposed  to  the  exaggerated  grains  of  a-SiC  which  can  grow  in 
sintered  8 microstructures.  ’ At  2055®C  (Figure  7e)  the  fine- 
grained matrix  has  disappeared  but  the  microstructure  is  still 
primarily  8-SiC  (approximately  5 percent  4H-SiC  as  determined  by 

4 

quantitative  X-ray  diffraction  ) . By  2105®C  (Figure  7F) , large 
portions  have  transformed  into  100-500  pm  "feather"  structures  of 
o-SiC.^  The  area  photographed  in  Figure  7f,  however,  is  a region 
between  "feathers"  which  remains  primarily  B-SiC. 

The  choice  of  1940®C  for  the  refiring  of  the  hot-pressed  bend 
specimens  was  based  on  both  this  microstructural  study  and  the 
refiring  data  included  in  Section  III  and  represents  a compromise 
between  a minimal  change  in  microstructure  (1940®C  resulted  in  a 
microstructure  between  Figures  33b  and  c)  and  the  less-than- 
optimum  crack  "blunting"  treatment  at  1940®C  as  illustrated  earlier 
in  Figure  2. 
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b.  1920“ 

Figure  33.  (Coni' d) 


a.  As  Hot  Pressed 


Figure  33. 


Optical  Micrographs  of  Polished  and  Etched 
Samples  of  Hot-Pressed  SiC  Annealed  for 
30  Minutes.  1250X,  dual  etch-molten  salt 
and  Murakami  reagent. 


d.  2015“ 

Figure  33.  (Cont'd) 


nir?' 


Figure  33.  (Cont'd) 


As  discussed  in  more  detail  in  Reference  3,  the  combined 

molten  salt  and  Murakami  etches  reveal  useful  detail  of  partially 

transformed  microstructures  that  cannot  be  resolved  by  either  etchant 

when  used  alone.  Transformation  of  3 to  o commonly  proceeds  in  sintered 

Sic  as  the  growth  of  thin,  tabular  grains  of  a-SiC  oriented  with  the 

basal  plane  of  the  crystallographic  structure  parallel  to  the  plane 

of  the  tabular  grains.  Careful  transmission  electron  microscopy  by 

25 

Shinozaki  and  Kinsman  first  identified  an  envelope  of  epitaxially 
connected  3-SiC  surrounding  each  tabular  grain  of  a-SiC  with  a (111) 
plane  of  the  cubic  structure  parallel  to  the  basal  plane  of  the  a-SiC 
and  with  a highly  planar  boundary  connecting  the  two  polytypes.  As 
reported  in  Reference  3,  the  molten  salt  etch  primarily  reveals  the 
irregular  external  boundary  where  the  envelope  of  the  B-SiC  contacts 
the  matrix  of  smaller  3-SiC  grains.  The  Murakami  etch,  on  the  other 
hand,  primarily  reveals  the  epitaxial  boundary  between  the  two  polytypes. 
The  combined  etch  of  partially  transformed  sintered  SiC  therefore  reveals 
both  the  planar  a- 3 boundaries  and  the  irregular  3-3  boundaries. 

Similar  features  can  be  resolved  in  the  microstructures  of  the 
hot-pressed  and  annealed  specimens  discussed  above  and  included  as 
Figures  33a-f.  Throughout  this  series  of  micrographs,  it  is  evident 
that  each  of  the  exaggerated  grains  of  "primarily  B-SiC"  contains 
a thin  layer  with  nearly  planar  boundaries.  These  are  particularly 
obvious  in  Figures  33c  and  d but  can  be  resolved  even  in  the  as-hot- 
pressed  microstructure  of  Figure  33a.  Selected  area  electron  diffraction 
and  X-ray  diffraction  have  identified  the  planar  layers  as  primarily 
4H-SiC  and  the  envelopes  as  3C-SiC  oriented  in  the  manner  discussed 
above.  The  shape  of  the  resulting  exaggerated  grains  in  fine-grained 
hot-pressed  SiC  is  closer  to  equiaxed  (thicker  layer  of  3-SiC  envelope) 
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grains,  the  rate  of  three-dimensional  grain  growth  is  probably  con- 
trolled by  the  rate  of  lateral  growth  of  the  a-SiC  lay.-r.  The  micro- 
structure of  the  hot-pressed  SiC  can  therefore  be  strongly  influenced 
by  small  amounts  of  0 to  a transformation.  The  total  amount  of 


i 


transformation  can  be  as  little  as  two  percent  and  result  in  a micro- 
structure such  as  Figure  33c  or  as  little  as  five  percent  and  ^esult  in 
a microstructure  such  as  Figure  33e. 
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X. 


INCREASING  FRACTURE  TOUGHNESS  THROUGH  TOE  ADDITION  OF  SECOND  PHASES 


Any  method  by  which  the  fracture  toughness  or  fracture  energy 

of  a ceramic  can  be  increased  without  increasing  the  inherent  flaw 

size  will  be  beneficial  with  respect  to  several  strength  related 

properties  including  the  fracture  strength  (through  the  Griffith 

relation)  and  resistance  to  impact  damage  (both  Hertzian  impact  and 

machining  damage) . There  has  been  a great  deal  of  work  reported  in 

the  literature  with  the  goal  of  increasing  the  fracture  toughness 

of  ceramics  through  the  addition  of  dispersions  of  second  phase 

particles.  This  area  has  been  reviewed  in  some  detail  recently  by 

26 

Donald  and  McMillan.  The  actual  role  of  second  phase  particles 

is  very  complex  and  depends  on  factors  such  as  thermal  expansion 

mismatch  (between  particles  and  matrix  phases),  elastic  modulus 

mismatch,  particle  size,  etc.  Depending  on  the  parameters  listed 

above,  the  fracture  toughness  is  generally  hoped  to  increase  through 

27 

interactions  of  the  crack  front  with  the  second  phase  particles, 

28 

with  stress  fields  in  the  matrix  surrounding  the  particles  or 
with  regions  of  microcracks  caused  by  the  particles.  The  actual 
interactions  may  cause  the  crack  to  be  arrested,  branched  or  deflected 
with  the  desired  net  effect  of  an  increase  in  the  energy  needed  for 
a unit  of  crack  propagation. 

During  the  previous  year,  several  preliminary  experiments  were 
directed  toward  the  study  of  particulate  second  phase  additions  in 
boron-doped  sintered  SiC.  The  eventual  goal  of  these  experiments 
is  a practical  means  of  increasing  the  fracture  toughness  of  sintered 
SiC  with  negligible  degradation  of  other  important  properties  such 
as  creep,  stress  rupture  and  oxidation.  In  order  to  accomplish 
this  goal,  several  criteria  must  be  satisfied  in  addition  to  those 
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which  directly  affect  the  fracture  toughness.  The  additional 
criteria  relate  to  the  compatibility  of  the  second  phase  particles 
with  the  matrix  SiC.  The  second  phase  additive: 

e must  not  inhibit  the  sinterability  of  boron-doped 
Sic  powder  through  either  chemical  or  physical 
mechansims, 

• must  exist  in  the  sintered  microstructure  as  isolated 
second  phase  particles  rather  than  a continuous  grain 
boundary  phase  (which  may  degrade  creep  and  stress  rup- 
ture properties) , 

• must  not  adversely  affect  the  protective  nature  of  the 
silica  layer  through  chemical  "fluxing",  spalling,  etc. 

Several  promising  second  phase  systems  were  identified  for 
this  preliminary  study  including  boron  carbide  and  transition  metal  ' 

diborides  with  volatile  oxides  such  as  M0B2  and  WB2.  The  physical 
and  chemical  properties  of  these  materials  combined  with  our  present 
understanding  of  the  sintering  process  in  boron-doped  SiC  suggested 
that  they  had  the  potential  of  meeting  the  compatibility  criteria 
listed  above.  The  additives  were  generally  in  the  form  of  1 to  3 pm 
particles  and  were  ball  milled  with  the  SiC  powder  for  4 to  6 hours 
to  provide  intimate  mixing.  The  resulting  microstructures  generally 
revealed  the  second  phase  particles  as  isolated  inclusions  with 
approximately  the  same  distribution  of  sizes  as  that  present  in  the 
starting  powder.  It  was  also  generally  observed  that  larger  amounts 
of  additive  enhanced  the  8 to  a polytype  transformation. 

To  date,  the  majority  of  work  on  second  phase  additives  in  SiC  f 

I 

has  been  aimed  at  studies  of  sinterability  using  approximately  0.5 
in.  diameter  pellets.  The  only  measurements  of  fracture  toughness 
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were  carried  out  on  two  groiq>s  of  double  torsion  specimens  containing 
approximately  0.5  and  5.0  w/o  of  boron  carbide  respectively.  The 
measurements  of  at  room  temperature  revealed  no  distinguishable 
difference  for  the  two  levels  of  additive  with  all  Kj^,  values  in 
the  vicinity  of  2.7  Kpsi  "VTii. 
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STRESS  RUPTURE  TESTING  OF  SILICON  CARBIDE 
AT  VERY  HIGH  TEMPERATURES 

R.J.  Charles  and  S.  Prochazka 
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INTRODUCTION 

As  ceramic  materials  continue  to  be  developed 
for  purposes  of  sustaining  mechanical  load  at  very 
high  temperatures,  the  problem  of  categorizing  the 
load  bearing  life-time  performance  characteristics 
of  these  materials  becomes  increasingly  difficult 
and  experimentally  time  consuming.  Two  major 
factors  contribute  to  such  evaluations.  On  the  one 
hand,  the  temperatures  involved  are  extreme  if 
rupture  at  nominal  stresses  is  to  be  achieved  in 
reasonable  times;  and  second,  the  volumes  of  ma- 
terial under  desired  tensile  stresses  are  likely  to 
be  small  if  the  normal  and  convenient  bend  test  ex- 
perimentation is  to  be  utilized. 

The  latter  factor  is  of  great  importance  because, 
if  the  results  are  to  be  used  for  components  of  prac- 
tical size,  a reasonably  accurate  extrapolation  of 
strength  capability  as  a function  of  stressed  volume 
must  be  undertaken.  Clearly,  the  accuracy  of  this 
extrapolation  is  improved  if  the  degree  of  extrapola- 
tion is  reduced.  Sedlacek  and  Halden^D  have  devel- 
oped a "pressurized  tube  technique"  for  the  strength 
testing  of  ceramics  at  room  temperature.  This  test 
offers  the  advantage  of  a large  stressed  volume  re- 
sulting from  the  nearly  uniform  tensile  stresses  gen- 
erated throughout  the  wall  of  the  pressurized  tube. 

In  fact,  with  tubes  of  reasonable  configuration,  the 
effective  stressed  volumes  can  be  increased  by  sev- 
eral orders  of  magnitude  compared  to  those  encoun- 
tered with  the  same  amount  of  specimen  material  in 
a bend  test  arrangement. 

The  object  of  the  present  work  is  to  report  a 
pressurized-tube,  stress-rupture  experimental  pro- 
cedure, using  sintered  silicon  carbide  material  at  a 
very  high  temperature,  and  to  correlate  the  results 
with  those  theoretically  expected. 

EXPERIMENTAL  PROCEDURE 

Closed-end  tubes  26  cm  long,  approximately 
1. 3 cm  OD,  and  1. 0 cm  ID  (in  the  final  fired  state) 
were  slip-cast  from  an  aqueous  dispersion  of  silicon 
carbide  stabilized  by  tetramethylammonium  hydroxide 
at  pH  10  and  1/2%  by  weight  of  solid  of  urea-formal- 
dehyde resol.  + Slip-casting  was  done  in  plaster-of- 
paris  molds,  and  the  castings  were  fired  at  2150  ®C 
in  1 atm  of  argon  to  a final  density  of  about  96%. 

The  closed  ends  of  the  tubes  were  hemispherical 
with  wall  thicknesses  of  about  1.5  mm. 

Pressurizing  of  a thin-walled,  closed-end  tube 
results  in  tangential  and  longitudinal  tensile  stresses 
of: 

♦URAC  180,  American  Cyanamide  Company. 


respectively,  where  P is  the  inside  hydrostatic  pres- 
sure, d is  the  diameter,  and  t is  the  wall  thickness. 

These  equations  are  strictly  correct  only  for 
thin-walled  pressure  vessels  where  the  inside  diam- 
eter is  essentially  equivalent  to  the  outside  diameter. 
The  ratio  of  wall  thicknesses  to  diameters  of  the  tubes 
in  the  present  tests  would  not  generally  be  classified 
as  those  for  a thin-walled  vessel,  but  for  the  sake  of 
simplicity  in  this  initial  study  such  an  assumption 
was  made.  The  tubes  were  circular  within  0.  2 mm. 
Therefore,  no  correction  for  ovality  was  necessary. 

The  open  end  of  the  tube  was  sealed  into  a pres- 
sure head  with  epoxy  resin  and  "proof-tested"  at  room 
temperature  by  pressurizing  to  15.  2 MPa  (2200  psi) 
nitrogen  pressure  which  resulted  in  hoop  stresses  of 
^approximately  51. 5 MPa  (7480  psi).  The  tube  was 
then  inserted  into  a graphite  resistance  laboratory 
furnace  of  the  type  described  by  St.  Pierre  and 
Curran.  (2)  Gas  cylinders  were  used  as  the  source 
of  nitrogen,  and  the  pressure  was  controlled  with  a 
pressure  regulator. 

The  temperature  was  measured  by  an  optical 
pyrometer  sighting  on  the  closed  end  of  the  tube 
through  a window  in  the  opposite  end  of  the  furnace. 

By  adjusting  the  position  of  the  tube  progressively 
farther  into  the  furnace,  it  was  possible  to  measure 
the  temperature  gradient  near  the  hot  zone  and  to 
position  the  tube  such  that  the  hipest  temperature 
was  about  1 cm  from  the  end.  In  this  position  there 
was  a temperature  drop  from  the  hottest  region  to 
the  closed  end  of  about  30  °C  over  a total  length  of 
approximately  1 cm.  The  temperature  as  measured 
by  the  pyrometer  was  corrected  for  losses  by  the 
furnace  window  and  a mirror.  An  accuracy  of  ±20  “C 
is  estimated  for  the  quoted  temperatures,  and  the 
region  of  highest  temperature  is  estimated  to  be 
2 cm  long. 

Numerous  safety  precautions  were  taken  to  avoid 
injury  at  the  time  of  specimen  failure.  Rather  than 
position  the  optical  pyrometer  directly  outside  the 
furnace,  a mirror  was  used  allowing  the  pyrometer 
and  the  operator  to  be  positioned  behind  a barrier. 
The  circumference  and  both  ends  of  the  furnace  were 
shielded  to  help  contain  any  fragments  of  hot  speci- 
men, etc.  which  might  be  released  if  the  furnace 
shell  should  fail.  Despite  such  precautions,  it  was 
estimated  that  the  graphite  felt  insulation  inside  the 
furnace  would  be  of  sufficient  thickness  to  absorb  the 
energy  and  contain  the  specimen  fragments,  thus 
preventing  damage  to  the  furnace  shell.  It  was  ob- 
served during  the  experiments  that  fragments  from 
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failed  specimens  were  contained  within  the  furnace, 
requiring  only  new  insulation  and  a new  heating  ele- 
ment to  repair  the  damage.  Nevertheless,  similar 
safety  precautions  employing  barriers  are  strongly 
recommended  in  the  use  of  such  a testing  procedure. 

Table  I lists  the  temperatures,  wall  stresses, 
and  test  times  for  tests  on  three  tubes  as  described 
above. 

TABLE  I 


Experimental  Tube  Bursting  Data 


Tube 

Temp 

Tangential 
Wall  Stress 

Time 

Test 

UK) 

(MPa/kpsi) 

(min) 

A 

2303 

35.2/5. 11 

90 

2353 

35.2/5. 11 

105 

2373 

35.2/5.11 

45 

2283 

46.4/6.  74 

30 

2303 

46.  4/6.  74 

15 

2353 

46.4/6.  74 

15 

2383 

46.4/6.  74 

23  (failure) 

B 

2433 

28.1/4.08 

15 

2433 

56.2/8.16 

0.  25  (failure) 

C 

2453 

32,0/4.  64 

12  (failure) 

THEORY  AND  ANALYSIS  OF  RESULTS 


Prochazka,  Johnson,  and  Giddings^^)  have  shown 
that  die-pressed  and  sintered  forms  of  silicon  car- 
bide, processed  for  high  average  strength,  exhibit 
Weibull  moduli,  m.  of  around  S to  6.  By  using  the 
two-parameter  Weibull  formulation,  we  may  esti- 
mate the  "quick  time"  or  low-temperature  bursting 
strength  of  that  section  of  the  tube  (>»  1 inch  long)  in 
which  stress  rupture  failure  would  be  expected  to 
occur  at  high  temperature  in  the  previously  described 
experiments.  The  two-parameter  Weibull  formula 
may  be  given  as  in  Eq.  (1 ),  where  P is  the  cumulative 
probability  of  failure  of  a given  structure;  om  is  the 
maximum  tensile  stress  in  the  structure;  KV  is  the 
"effective  volume"  under  tension;  and  oo  is  a nor- 
malizing constant. 

tn  tn  ^ I'p  « m -In  Om  + fn  (KV/ctq*")  . (1) 

For  three-point  bend  K « 0.  5/(m  + 1)*  and  pure  tension 
K ■■  1;  thus,  for  a constant  probability  of  failure  of 
two  structures  where  one  is  a standard  three-point 
bend  specimen  and  the  other  is  a volume  of  material 
under  uniform  tension,  one  has 

<Jt/0B  * <0-5(VB/Vx)/(m+l)*)*/'"  (2) 


estimation  we  use  a value  of  5 consistent  with  values 
appropriate  for  die-pressed  materials.  As  reported 
previously<3)  the  bend  bar  shape  utilized  in  testing  of 
these  materials  is  0.  4 cm  by  0. 4 cm  by  4.  0 cm 
(0. 16  inch  by  0. 16  inch  by  1.6  inch).  Thus,  for  the 
tubes  and  bend  bar  geometries  the  ratio  aj/oB  cal- 
culates to  be  about  0.35  — and  for  an  MOR  of  4.13 
MPa  (60  kpsi),  appropriate  to  the  slip  cast  material 
as  shown  by  Prochazka,  (4)  the  expected  "quick  time" 
or  low-temperature  rupture  stress  for  the  slip  cast 
tube  end  (1  inch  long)  would  be  about  1.43  MPA 
(20.  7 kpsi). 


An  estimate  of  the  stress  rupture  character- 
istics of  dense,  sintered  silicon  carbide*  at  very 
high  temperatures  has  been  recently  made(^)  by  a 
theoretical  extrapolation  of  1 600  °C  rupture  data  given 
for  this  material  by  Trantina  and  Jo)ui8on(^)  and  by 
Prochazka.  C?)  The  theory(8)  for  the  extrapolation 
was  based  on  the  concept  of  the  growth  in  the  material 
of  preexisting  flaws  to  critical  sizes  which  would  re- 
sult in  catastrophic  failure  under  the  applied  load. 
Crack  growth  was  related  to  diffusional  processes  in 
the  material  such  that  the  formulation  for  predicting 
stress  rupture  behavior  involved  a parameter  which 
could  be  described  as  a "temperature  modified  failure 
time"  (i.  e. , the  time  integral  of  the  exponential  of  an 
activation  energy  of  diffusion  divided  by  temperature 
and  the  Boltzmann  constant).  Such  a parameter  is 
similar  to  the  Orr-Sherby-Dorn  (OSD)  parameter,  (9) 
applicable  in  the  stress  rupture  analysis  of  metals, 
and  is  referred  to  in  the  present  work  by  the  same 
symbol,  Qj.#  utilized  earlier  by  these  authors. 


Utilizing  this  parameter,  the  result  of  the  stress 
rupture  theory  indicated  above(8)  may  be  given  in  the 
following  mathematical  form: 


Z(N,  R)  = / 


,tr  exp  (-AHc/kT)dt  6r 


00 


= 7^  (3) 


In  Eq.  (3),  tj.  is  the  rupture  time,  R is  the  ratio 
Oa/omax  where  Oa  is  constant  stress  applied  to  the 
sample,  and  Omax  is  an  upper  limit  of  stress,  which 
when  applied,  results  in  immediate  failure.  The  latter 
stress  is  thus  independent  of  diffusional  processes  and 
reflects  the  strength  limiting  flaws  in  the  as-processed 
material.  For  a given  material  or  sample,  this  stress 
may  most  conveniently  be  determined  by  tests  at  suf- 
ficiently low  temperatures  that  diffusional  effects 
would  be  negligible.  Also  in  Eq.  (3),  kT  has  its  usual 
significance;  AHc  is  an  activation  energy  for  creep  or 
the  diffusion  of  those  species  responsible  for  creep; 

6o  and  N are  constants  for  a given  material,  which 
may  be  estimated  theoretically  but  are  more  appro- 
priately determined  experimentally;  and  Z(N,  R)  is  a 
stress  function  which  results  from  the  theory  to  be  of 
the  following  form. 


where  the  subscripts  T and  B refer,  respectively,  to 
the  structure  under  uniform  tension  and  the  three- 
point  bend  specimen.  For  the  slip-cast  sintered  SiC 
mrterial  the  appropriate  Weibull  modulus,  m,  has 
not  been  determined.  Therefore,  for  purposes  of 


Z(N.  R) 


e-N(N’  + 3n4  + 6N  + 6) 


(4) 


*Fractional  % boron  additives  utilized  as  sintering  aids. 
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It  may  be  noted  in  Eq. 
is  constant  then 


(3)  that  if  the  temperature 


Gr  = tr  exp  (-iHc/kT) 


(5) 


and,  as  such,  is  identical  with  the  definition  of  this 
parameter  as  originally  given  by  OSD.  However,  if 
the  temperature  varies  with  time  during  the  test  then 
6r  is  a summed  or  integral  value  as  indicated  in 
Eq.  (3).  It  will  be  necessary  to  use  the  more  gen- 
eral definition  in  the  following  analysis  of  tube  burst- 
ing data  because  temperature  changes  in  fact,  did, 
occur  during  the  course  of  a single  experiment. 

It  should  further  be  noted  that  Eq.  (3)  applies 
to  a simple  rupture  experiment  (a  single  constant 
stress)  and  involves  an  approximation  arising  from 
a disregard  of  a limit  of  intergration  for  terms  on 
the  left-hand  side  of  the  equation.  (8)  As  such,  Eq.  (3) 
introduces  an  error  which  becomes  appreciable  as 
R approaches  the  numerical  value  of  one.  The  pres- 
ent work  will  be  concerned  with  varied  and  progres- 
sive changes  in  applied  stress  during  the  course  of 
single  rupture  experiment.  Thus,  it  will  be  neces- 
sary to  account  for  the  effects  of  prior  stress  (and 
temperature)  histories  on  the  growth  of  cracks  and 
the  progressive  weakening  of  a given  test  sample. 

In  such  cases,  neglect  of  the  integration  limit,  as  is 
appropriate  for  most,  simple  rupture  experiments, 
is  not  proper.  With  both  limits,  Eq.  (3)  reverts  to 
the  following  equation  which  can  properly  account 
for  interim  states  in  the  process  toward  rupture, 

J**"  exp  (-AHc/kT)  dt 

Z(N,  Ri)  - Z(N,  Rq)  =-2 


_0_ 

00 


(6) 


TABLE  n(a) 


Calculation  of  9r 

for  Tube  A 

Temp 

Ri 

9 

(°K) 

(hr  X 10 

2303 

0.247 

0.40 

2353 

0.247 

4.02 

2373 

0.  247 

14.85 

2283 

0.247 

9.96 

re  * 2,  92  X 10-*’  hr 


re 


Z(N,  Rq)  •=  Z(N,  Rj)  - 


Z(39,  0.247) 
5.  6 X 10” 


2.  92  X 10'** 
5. 12  X 10"* 


Ri  is  the  initial  value  of  Oa/Omax  Rq  in- 

terim value.  It  may  be  well  to  point  out  that  the  dif- 
ference between  Ri  and  Rq  is  the  result  of  the  change 
in  Omax>  ih®  quick  time  strength  of  the  sample,  under 
the  constant  applied  stress,  oa*  This  change  results 
from  the  growth  of  a strength  limiting  flaw  from  its 
initial  size  to  some  interim  size  with  a consequent 
reduction  in  instantaneous  strength  of  the  sample. 

Equation  (6)  is  the  appropriate  equation  for  ana- 
lyzing the  data  in  Table  I to  determine  if  the  prelim- 
inary and  lower  stress  applications  contributed  any 
significant  weakening  of  tlie  specimen  before  the 
higher  and  final  stress  was  applied. 

From  other  work^S)  a value  of  AHc  « 250  kcal/mol 
and  a log  do  ■ -35.  33  have  been  indicated  as  appro- 
priate for  dense,  sintered  silicon  carbide  and,  since 
in  the  previous  discussion  a value  of  0.  35  was  cal- 
culated for  the  ratio  Of/oB,  the  expected  quick  time 
strength,  Oinax>  ^h®  slip-cast  burst  tubes  would 
be  0.  35  (41 3.  3)  « 143  MPa  (20.  7 kpsi).  Table  H lists 
a two- stage  calculation  to  determine  the  appropriate 
8r  value  for  tube  A. 


Rq  = 0.  249 

We  conclude  from  Table  11(a)  that  the  preliminary 
stress  and  thermal  history  for  tube  A results  in  a 
slight  weakening  of  the  tube  amounting  to  about  1%  of 
its  original  strength  capability  [i.  e. , (.  249  - . 247)/ 

. 247  ~ 0.  01].  Thus  the  apparent  stress  ratio  R(0.  326) 
given  in  Table  I for  the  remainder  of  the  experiment 
should  be  increasecf  by  this  small  factor  [ i.  e. , 1.01 
(.  326)  ~ 0.  329],  but  it  is  clear  that  such  a correction 
is  small  relative  to  other  uncertainties  in  the  data. 

The  above  calculation  does  indicate,  however,  that 
the  lifetime  of  the  sample  is  governed  almost  totally 
by  crack  growth  under  the  final  stress  at  which  rupture 
was  observed.  Table  11(b)  gives  the  remaining  calcu- 
lation for  for  tube  A under  these  conditions. 

TABLE  n(b) 

Calculation  of  0^.  for  Tube  A 


Temp 

(°K) 

2283 

2303 

2353 

2383 


Ri 

(oa/cmax)i 

.329 
. 329 
.329 
.329 


(hr  X 10**) 

0.83 
0.  67 
2.12 
6.  35 


re  = 9j.  = 9.  97xl0-*‘  hr 


Similar  calculations,  accounting  for  the  prior 
stress  history  of  tube  B,  indicate  that  virtually  no 
weakening  of  this  tube  occurred  as  a result  of  the  pre- 
liminary loading  and  that  rupture  was  almost  totally 
determined  by  the  final  a.2d  higher  stress.  Table  III 
therefore  summarizes  the  summed  or  integral  OSD 
rupture  parameters  for  all  three  tubes. 

Figure  1 plots  the  data  (solid  circles)  listed  in 
Table  in  on  the  estimated  stress  rupture  curve  for 
sintered  silicon  carbide  which  was  prepared  earlier 
(Ref.  5)  from  rupture  data^®»  obtained  at  1600  °C 
(open  circles)  and  the  equation  referred  to  in  this  work 
as  Eq.  (3).  It  is  evident  that  the  stress  rupture  results 
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Fig.  1 Estimated  stress  rupture  curve  for  dense, 
pure-phase  sintered  silicon  carbide  according  to 
Eq.  '3)  and  the  physical  parameters  listed  in  the 
figure.  Open-circle  points  calculated  from  exper- 
imental bend  bar  data  at  1600  closed-circle 

points  calculated  from  experimental  data  of  this 
work. 


TABLE  m 

Summary  of  6p  Values 


Tube 

Test 

R 

(Oa/Omax) 

fir 

(hr  X 10“) 

6r* 

A 

0.  329 

9.9 

12.  33 

B 

0.  394 

0.203 

10.  64 

C 

♦Log  9o 

0.224 

“35.33(5) 

14.8 

12.50 

of  this  work,  obtained  at  low  stresses  and  very  high 
temperatures,  are  reasonably  consistent  with  the 
originally  estimated  characteristics  of  this  material. 

It  should  be  noted,  however,  that  the  approach 
described  herein  involves  several  uncertainties  and 
estimates  related  to  the  data.  In  particular  some 
concern  exists  that  the  two-parameter  Weibull  for- 
mulation. which  is  used  in  the  present  work  to  cor- 
relate sample  strength  with  stressed  volume,  may 
overestimate  the  degradation  in  strength  because  of 
increased  volume  under  stress.  (3)  On  the  other  hand, 
the  bursting  behavior  of  the  tubes  was  correlated  to 
only  the  tangential  stress  in  the  tube  walls.  The  true 
stress  state  is  biaxial  and  the  neglect  of  the  action 
of  flaws  brought  into  unfavorable  orientation  with 


♦Earlier  analyses  of  crack  dynamics  governing  de- 
layed failure  identify  this  same  time  integral  as  the 
determining  factor  in  rupture  life  in  systems  involv- 
ing thermal  activation  (e.  g. . Ref.  10). 


respect  to  the  smaller,  longitudinal  stress  would 
compensate,  to  some  degree,  an  overestimation 
of  strength  degradation  resulting  from  increased 
volume  under  stress. 

The  present  ainalysis  suggests  a further  aspect 
of  stress  rupture  testing  that  may  have  general  ap- 
plicability to  materials  in  which  diffusional  processes 
play  a dominant  role  in  time- dependent  failure.  For 
rupture  under  a given  stress,  the  parameter  govern- 
ing lifetime  is.  according  to  Eq.  (3).  the  time  in- 
tegral of  the  exponential  of  -AH/kT.  ♦ For  a constant 
rate  of  temperature  increase.  6.  the  relation  between 
the  integral  OSD  parameter.  6r.  and  the  rate  of  tem- 
perature increase  would  be  simply. 

T 

9r  ='|-[Texp(-AH_/kT)  + Ei(-aHc/kT]_*'  (7) 

B 10 

where  the  evaluation  limits  are  the  initial  tempera- 
ture. To.  and  the  higher  rupture  temperature.  Tj.. 
and  the  symbol  Ei  refers  to  tabulated  exponential 
integrals.  For  nonlinear  temperature  schedules, 
similar  integrals  may  be  evaluated.  It  is  evident 
that  data  necessary  to  describe  the  stress  rupture 
characteristics  of  a material  could  be  obtained  in  a 
dynamic  fashion  in  which  the  experimental  procedure 
would  be  to  raise  the  temperature  in  a controlled 
manner  and  note  the  final  temperature  at  which  rup- 
ture occurred.  This  procedure  would  introduce  a 
degree  of  control  in  developing  test  data  that  is  absent 
in  the  normal  procedure  where  stress  and  tempera- 
ture are  fixed,  and  there  is  no  assurance  that  failure 
can  be  achieved  within  a practical  time  interval. 

In  many  respects  the  above  procedure,  involving 
a variation  in  temperature,  is  complementary  to  a 
procedure  suggested  earlier(ll)  to  assure  that  failure 
occurred  within  practical  times.  This  latter  procedure 
involved  a controlled  increase  of  applied  stress  to  the 
rupture  point.  For  the  results  to  be  useful,  however, 
for  generalizing  on  stress  rupture,  it  was  necessary 
to  assume  a relation  between  crack  growth  velocity 
and  applied  stress.  Such  relations  are  exceedingly 
difficult  to  define  with  high  accuracy.  In  the  present 
scheme,  involving  temperature  variations,  it  should 
be  possible  to  determine  the  unknown  factor  in  Eq.  (7). 
aHc.  to  high  accuracy  with  only  a few  controlled 
measurements.  Once  the  factor.  AHc<  is  established 
then  experiments  at  different  stress  levels  could 
rapidly  develop  the  necessary  stress  dependencies 
of  failure  time. 

SUMMARY  AND  CONCLUSIONS 

A bursting  tube  testing  procedure,  described  by 
Sedlacek  and  Halden.  ) has  been  evaluated  as  a means 
of  obtaining  high-temperature  stress  rupture  data  on 
pure-phase.  sintered  silicon  carbide.  This  procedure 
was  selected  since  it  offered  a simple  method  for  high- 
temperature  experimentation  and  the  possibility  of  in- 
creasing the  effective  volume  under  stress  by  a factor 
of  about  two  orders  of  magnitude  greater  than  that  ob- 
tained in  small  bend  bar  testing.  Analysis  of  the  re- 
sults indicates  that: 


i 

1 
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1 . The  bursting  tube  technique  is  a useful  and 
convenient  procedure  for  the  extension  of  high- 
temperature  mechanical  properties  measure- 
ments to  specimens  of  relatively  large  stressed 
volume. 

2.  In  the  range  of  2000  °C,  dense,  pure-phase 
sintered  silicon  carbide  exhibits  practical 
strengths  and  tnipture  life  times  in  inert 
atmospheres. 

3.  Parameters  have  been  defined  for  a stress 
rupture  function  which,  when  combined  with  a 
simple  form  of  Weibull  statistics,  permits 
reasonable  estimates  of  the  stress  rupture 
performance  of  this  material  at  high  tempera- 
ture in  varied  sample  configurations. 

4.  Stress  rupture  testing  procedures  which  employ 
controlled  temperature  variations  are  analyzable 
and  may  have  general  applicability. 
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mCH-TEMPERATURE  STRESS  RUPTURE 
OF  POLYCRYSTALLINE  CERAMICS 


R.  J.  Charles 


INTRODUCTION 

An  analysis  of  the  high- temperature  stress  rup- 
ture process  in  polycrystailine  materials  has  recently 
been  presented(l)  in  which  the  major  factors  gov- 
erning time- dependent  rupture  were  identified  as  the 
diffusivities  of  the  atomic  species  responsible  for 
diffusional  creep  and  the  dependence  of  these  diffusiv- 
ities on  mechanical  pressures.  Creep  at  hi^  temper- 
atures was  assumed  to  result  from  the  processes  de- 
scribed by  Nabarro'^^  and  Herring^^'  for  the  relax- 
ation of  stress  by  the  diffusion  of  atoms  or  vacancies 
from  one  pressure  state  to  another.  A fracture  me- 
chanics approach  was  utilized  to  describe  these  pres- 
sure states,  such  that  flaws  within  grain  boundaries 
would  propagate  under  an  applied  load  and  ultimately 
result  in  rupture.  Thus,  in  principle,  the  stress  rup- 
ture performance  of  a polycrystailine  material  could 
be  described  wholly  in  terms  of  diffusional  parameters 
governing  creep,  the  pressure  and  temperature  depen- 
dencies of  these  parameters,  and  time-independent 
physical  properties  of  the  material. 

While  the  above  analysis  was  primarily  under- 
taken to  provide  information  which  would  assist  in 
categorizing  the  high-temperature  mechanical  prop- 
erties of  a number  of  recently  developed  structural 
ceramics,  it  was  foiand  that  little  data  existed  in  the 
literature  on  stress  rupture  of  ceramics  to  which  the 
above  analysis  covild  be  applied  for  purposes  of  veri- 
fication. An  abundance  of  stress  rupture  data  exists 
for  metals  and  alloys.  Application  of  the  analytical 
procedures  developed  in  the  preceding  work  to  a num- 
ber of  pure  metal  and  alloy  systems  indicated  that  the 
stress  rupture  characteristics  of  these  materials 
could  be  accoimted  for  with  good  accuracy  over  a large 
range  of  times,  temperatures,  and  applied  stresses. 
Similar  descriptions  should  be  possible  for  the  stress 
rupture  characteristics  of  polycrystailine  ceramics. 

The  objectives  of  the  present  work  are  to  apply 
the  above  procedures  to  the  limited  stress  rupture 
data  available  on  ceramics  and,  more  particularly, 
to  use  these  procedures  to  provide  estimates  of  the 
potential  service  life  times  of  recently  developed  sil- 
icon nitride  and  silicon  carbide  structural  ceramics. 

BACKGROUND  THEORY 

The  stress  rupture  relationship,  referred  to  in 
the  previous  section,  is  given  by  the  following  equa- 
tion; 

e-RN  [nS/r+3n»/R»+6N/R’+6/R*1 

«8  Dg  (Ban,ax>*  tf  (1) 


In  the  above  equation,  R is  a dimensionless  ratio, 
Oa/Omax*  where  Oa  is  the  stress  applied  to  a sample 
and  ffjnax  upper  limit  of  applied  stress  which 

results  in  instantaneous  rupture  without  any  time- 
delayed  effects.  This  latter  stress,  is  gen- 

erally to  be  obtained  by  tests  at  sufficiently  low  tem- 
peratures that  diffusional  processes  would  have  neg- 
ligible effects  on  the  failure  process  and  thus  this 
stress  reflects  the  characteristics  of  the  strength 
limiting  flaws  existing  in  the  material  prior  to  test. 
The  factor  N is  a further  dimensionless  quantity  given 
by  the  product  where  Ujh  is  the  cohesive  or  theo- 
retical strength  of  a grain  boundary  and  B is  a constant 
determined  by  the  pressure  coefficient  of  diffusion  of 
the  mobile  species  in  the  material.  As  has  been  shown 
(Ref.  1),  e may  be  estimated  for  uniaxial  stress  situa- 
tions from  the  activation  volume,  v*,  of  the  diffusing 
species  according  to  the  following  equation, 

8«»v*/3kT  (2) 

where  kT  has  its  usual  significance. 

Further,  in  Eq.  (1 ),  Dg  is  the  diffusion  coeffi- 
cient of  mobile  species  and  may  be  given  in  its  usual 
form, 

Dg  = D^  exp(-AHj,/kT)  (3) 

where  A H^,  is  the  activation  energy  for  diffusion  and 
Djj  is  the  preexponential  factor.  The  final  three  pa- 
rameters Y,  P,  and  tf,  are,  respectively,  the  surface 
energy  of  the  material,  the  radius  of  curvature  of  an 
advancing  flaw  tip  (estimated  to  be  of  atomic  dimen- 
sions), and  the  time  to  failure  after  application  of  the 
applied  stress. 

It  is  evident  that  many  of  the  parameters  given 
in  Eq.  (1)  will  be  difficult  to  determine,  even  for  es- 
timates, and  indeed  the  equation  involves  a number  of 
simplifying  approximations.  Consequently,  the  best 
utility  of  Eq.  (1)  may  reside  in  its  use  to  rationalize 
existing  stress  rupture  data  or  to  extrapolate  such 
data  outside  of  the  range  of  tests. 

Recognizing  that  Eq.  (1)  embodies  the  well-known 
Orr-Sherby-Dorn  (OSD)(4)  stress  rupture  parameter, 
6j.  = tf  exp(-AHj,/kT),  where  AHj,  is  the  activation  en- 
ergy for  creep,  it  is  useful  to  collect  remaining  terms 
specific  to  a given  material  into  an  analogous  parame- 
ter, 00.  which  has  units  of  time  similar  to  0j..  Thus, 
Eq.  (1)  may  take  the  form, 

e-'^(R-l\NVR+3N’/R*+6N/R*+6/R*)_  ^ _ ®r 
(N’+3N®-i«N+6)  00 
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where  both  sides  of  the  equation  represent  dimension- 
less ratios.  If  the  activation  energy  for  creep.  AHg, 
is  known  for  a polycrystalline  ceramic  then,  in  prin- 
ciple, two  experimental  values  of  failure  times  tj,  and 
their  respective  applied  stress  ratios,  R = ^a/Omax  • 


are  sufficient  to  define  6,  and  N in  the  above  equation. 


Once  these  parameters  are  defined,  failure  times  un- 
der stress  mpture  may  be  estimated  for  all  ranges  of 
stresses  and  temperatures,  and  the  slope  parameter, 
N,  may  be  compared  with  that  theoretically  expected. 


The  use  of  Eq.  (4)  presupposes  that  an  estimate 
or  test  value  of  the  instantaneous  strength  (or  low- 
temperature  strength),  ^max'  available.  If  such  is 
not  the  case,  then  Eq.  (4)  may  be  cast  in  a further 
form  which,  however,  when  utilized,  provides  no  cor- 
roboratory information  concerning  the  expected  slope 
parameter  value,  N,  for  a given  system.  This  form 
is  given  in  Eq.  (5)  wherein  one  sets  V = 
thus 

e 


e°aV  (V^/a^  + 3V*/ag2  + 6V/Og^+  6/0^4)  = Z'-— p (5) 


where  is  a new  constant  for  a given  system. 


APPLICATIONS 


AlO, 
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Fig.  1 Stress  rupture  of  alumina  rods 
[after  Walles(5)]. 


a number  of  temperatures,  times,  and  applied  stresses. 
For  this  material,  low-temperature  results  were  also 
given  which  allow  estimates  of  the  instantaneous 
strength  of  the  fabricated  rods.  Consequently,  Eq.  (4) 
may  be  applied  to  these  data.  In  addition,  Rosolowski 
and  Greskovich^®)  have  summarized  a large  amoimt 
of  data  on  the  sintering  of  alumina,  and  conclude  that 
the  activation  energy  of  the  sintering  process  for  this 
material  is  most  characteristically  given  by  a value 
of  about  115  kcal/mole.  We  utilize  this  value  as  the 
appropriate  value  for  AHg  which  would  govern  the 
diffusion  mechanisms  for  both  sintering  and  creep  in 
polycrystalline  alumina.  Figure  2 gives  the  results 
of  applying  Eq.  (4)  to  Walles'  data  and  using  a value 
of  AHc  = 115  kcal/mole.  The  points  are  determined 
from  the  experimental  values  of  6j.,  utilizing  an  ar- 
bitrary constant  0 . whereas  the  solid  line  throu^ 
the  points  corres^nds  to  the  theoretical  function  on 
the  left-hand  side  of  Eq.  (4)  using  a slope  parameter 
(N)  value  of  16.  We  may  estimate  whether  or  not 
such  an  N value  is  reasonable  using  the  relation 
N ■ V*  oth/kT  and  setting  v*  as  about  5 cc/mole^Oth 
as  about  0.  2 of  Young's  Modulus  (3.8  x lO®  MN/m*) 
and  a temperature  of  about  1000  "C,  A calculated 
value  for  N of  about  12  results,  which  is  roughly  in 
accord  with  that  observed. 


Walles^®^  has  reported  several  series  of  stress  rup- 
ture results  on  fibers  and  rods  of  sintered  a.umina.  Fig- 
ure 1 reproduces  Walles'  results  for  alumina  for 


Fig.  2 Experimental  stress  rupture  values  of  Walles 
IRef.  (5)]  (points)  for  sintered  alumina  rods  compared 
with  the  stress  rupture  curve  (solid  line)  predicted 
by  Eq.  (4),  utilizing  the  physical  parameters  listed 
in  the  figure. 


As  indicated  in  Fig.  2 there  is  an  apparent  devia- 
tion of  experiment  from  theory  for  R values  less  than 
about  0.  2.  This  deviation  may  not  be  real  since  at 
such  low  R values  (low  applied  stresses)  the  crack 
sizes,  just  prior  to  rupture,  can  be  calculated  from 
the  Griffith  theory  to  be  a significant  fraction  of  the 
section  sizes  of  the  relatively  small  diameter  fibers 
used  in  the  experiments.  Consequently,  the  true 
stresses  on  the  unruptured  portions  of  the  fibers,  just 
prior  to  failure,  wotdd  be  significantly  larger  than  the 
apparent  stresses,  and  the  6r  points  would  therefore 
fall  below  the  expected  solid  line. 
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Walles  does  not  report  low-temperature  strength 
values  of  the  silicon  carbide  fibers  used  in  his  stress 
ruptiure  experiments  and  thus,  since  a value  of  the 
instantaneous  strength  is  not  available,  Eq.  (5)  is 
utilized  to  analyze  his  results.  Figure  3 reproduces 
the  original  test  data  and  Fig.  4 shows  the  results  of 
applying  Eq.  (5)  to  the  data.  In  this  case,  curves 
corresponding  to  the  function  represented  by  the 
left-hand  side  of  Eq.  (S)  were  first  passed  through 
the  experimental  data  for  each  temperature  at  which 
experiments  were  reported.  This  procedure  served 
to  identify  a value  of  ftHc  of  130  kcal/mole  as  best 
fitting  the  overall  data.  The  results  for  all  tempera- 
tures were  then  combined  to  result  in  the  representa- 
tion given  in  Fig.  4 in  which  the  solid  line  is  the  final 
curve  corresponding  to  the  Z'  function  and  the  points 
are  derived  from  experimental  8r  values.  It  is  evi- 
dent that  the  curvature  in  the  experimental  data  can 
be  well  represented  by  the  theoretical  stress  rupture 
function.  However,  it  is  of  some  interest  that  the  ap- 
propriate AHc  value  for  these  silicon  carbide  fibers 
is  as  low  as  130  kcal/mole.  The  data  of  Ghoshtagore 
and  CobleC^)  would  suggest  a substantially  higher  value 
(250  to  300  kcal/mole)  for  the  activation  energy  for 
the  species  (carbon)  controlling  diffusive  creep  in  sil- 
icon carbide.  Further,  the  results  are  of  significant 
interest,  in  that  even  at  short  times  there  exists  almost 
a threefold  reduction  in  strength  of  the  fibers  when 
tested  at  800  *C  and  tested  at  the  still  relatively  low 
temperature  of  1100  "C. 


with  the  stress  rupture  curve  (solid  line)  predicted 
by  Eq.  (5)  utilizing  the  physical  parameters  listed 
in  the  figure . 


tungsten  is  about  130  kcal/mole,  which  is  equivalent 
to  that  for  the  stress  rupture  process  exhibited  by 
Walles'  fibers.  It  is  possible,  therefore,  that  both 
the  high-temperature  strength  and  creep  rupture  prop- 
erties of  WaUes'  fibers  are  not  determined  by  SiC  but 
by  a substrate  material  such  as  tungsten.  If  this 
could  be  substantiated,  it  would  indicate  that  the  stress 
rupture  analysis  procedure  utilized  in  this  work  has 
value  not  only  for  extrapolatory  uses  but  also  for  eval- 
uating the  behaviors  of  complex,  composite  systems. 


Fig.  3 Stress  rupture  of  silicon  carbide  fibers 
[after  Walles(5)]. 

Silicon  carbide  fibers  are  generally  prepared  by 
chemical  vapor  deposition  processes  and,  as  such, 
are  usually  of  high  purity.  As  will  be  shown,  the 
dense,  pure-phase  silicon  carbide  shows  high  strength 
in  bulk  form  at  temperatures  of  1600  *C  or  higher. 
Consequently,  Walles*  fibers  of  SIC  do  not  show  the 
high-temperature  stability  normally  expected  for  a 
high-purity  SiC.  Often,  CVD  fibers  are  grown  on  a 
refractory  metal  substrate,  e.g..  tungsten.  We  note 
that  the  activation  energy  for  self-diffusion  in 


Sintered  SiC 

Trantina  and  Johnson(8)  and  ProchazkadD  have 
reported  limited  stress  rupture  data  on  a form  of 
dense,  pure-phase  sintered  SiC  described  earlier  by 
Prochazka^^^  The  information  consisted  of  both 
static  fatigue  (l.e.,  constant  stress)  and  constant 
stress  rate  data.  The  former  authors  report  their 
their  data  over  a limited  raiwe  according  to  a 
semi -empirical  relationshipdO)  derived  for  brittle 
matet  ials  on  the  basis  of  an  early  form  of  fracture 
mechanics  theory  and  which  may  be  glve.n  as  follows: 

In  R = In  {oj  Omax)  * I**  K 


S 


(6) 


where  K and  n are  constants  specific  to  the  material. 
Charles^l)  has  shown  that  the  slope  parameter,  N, 
utilized  in  the  present  work  bears  the  following  rela- 
tionship to  the  constant,  n,  in  Eq.  (6)  for  relatively 
large  values  of  R: 

1 (7) 

alntf  n NR 

Thus,  since  the  data  of  Trantina  and  Johnsonwere 
obtained  for  R values  centering  around  about  0.  84-- 
and  their  value  for  n was  given  as  about  33 — the  ap- 
propriate N value  for  the  present  analysis  would  be 
about  33/0.  84~39.  Prochazka's  creep  rate  studiesdD 
on  this  material  indicate  a temperature  dependence  of 
near  steady-state  creep  consistent  with  an  activation 
energy  of  about  210  to  250  kcal/mole.  Such  a value 
is  near  the  lower  end  of  values  determined  by  Coble 
and  GhoshtagoreC^)  for  carbon  diffusion  in  single- 
crystal SiC  (250  to  300  kcal/mole).  Using  a value 
for  ^Hc  of  250  kcal/mole  and  N = 39,  Eq.  (4)  may  be 
applied  to  the  combined  stress  rupture  data  (Refs.  8 
and  11)  to  obtain  an  estimation  of  the  stress  rupture 
performance  of  dense,  sintered  SiC  over  extended 
ranges  of  temperature  and  stress.  Figure  5 shows 
the  experimental  data  plotted  as  reduced  Orr-Sherby- 
Dorn  parameters  (6r/0o  where  the  appropriate 
logioSo  value  for  this  material  = - 35.33)  through 
which  is  passed  a curve  (solid  line)  corresponding  to 
the  Z stress  rupture  function  for  N • 39.  Using  this 
curve  and  a value  of  AHq  = 250  kcal/mole,  rupture 
times  for  sintered  silicon  carbide  may  be  estimated 
for  any  combination  of  stress  or  temperature. 
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Fig.  5 Estimated  stress  rupture  curve  for  sintered 
silicon  carbide  defined  by  experimental  data  (points) 
and  Eq.  (4)  utilizing  the  physical  parameters  listed 
in  the  figure. 

In  general,  designers  find  stress  rupture  infor- 
mation to  be  most  conveniently  represented  in  a plot 


of  applied  stress  vs  temperature  for  various  life  times 
(e.g. , 100,  1000,  10, 000  hours,  etc.)  Figures  6(a) 
and  (b)  are  such  representations  calculated  directly 
from  Eq.  (4)  where  logioBo  “ ' 35.33  and  AHg  ■ 250 
kcal/mole.  It  is  to  be  noted  that  the  Omax  value  uti- 
lized in  these  calculations  was  340  MN/m*  (50  KPSI), 
which  represents  a value  shown  to  be  appropriate  by 
Johnson  and  Trantina  for  sintered  SiC  in  pure  tension 
in  air.  The  pure  tension  amax  value  is  an  appropriate 
value  for  the  direct  comparison  of  stress  rupture  per- 
formance of  ceramics  vs  alloys.  Such  a comparison 
IS  shown  in  Fig.  6 wherein  the  stress  rupture  data  of 
Betteridge(12)  and  Goldhoff,  ^^3)  as  anadyzed  by  Charles 
(Ref.1)  for  Nimonic  80 A superalloy,  is  included.  Also 
included  in  these  figures  are  estimated  stress  rup- 
ture data  on  hot- pressed  forms  of  SISN4  and  SiC, 

These  latter  data  are  obtained  from  the  original  work 
of  Lange^^^^  and  were  calculated  using  the  same  pro- 
cedures as  those  outlined  for  the  analysis  of  the  data 
of  Trantina  and  Johnson(8)  on  sintered  SiC.  In  the 
CEtse  of  the  hot-pressed  materials,  Lange  reported 
constant  stress  rate  data  at  several  temperatures, 
thus  allowing  a direct  evaluation  of  the  appropriate 
stress  rupture  temperature  dependence  parameters, 
AHc.  Finally,  in  dashed  lines,  an  estimated  curve 
for  an  advanced  superalloy,  currently  under  develop- 
ment iy/y'-i),  (15)  Is  also  included  for  comparison 
purposes. 


Fig.  6(a)  One-hour  stress  rupture  curves  for  sev- 
eral high- temperature  structural  ceramics  and  for 
an  early  (Nlomic  80A)  and  an  advanced  (DSy/y'-5) 
superalloy  as  determined  by  Eq.  (4)  and  available 
experimental  data. 

Several  features  of  these  comparisons  are  of  in- 
terest. For  example,  superalloys  are  superior  in 
bulk  strength  at  low  temperatures  and,  in  fact, 
strengths  of  1000  MN/m*  (150  KPSI)  for  these  mate- 
rials at  low  temperatures  are  common.  However, 
there  is  a crossover  at  high  temperatures  where  the 
stress  rupture  performances  of  the  ceramics  exceed 
the  alloys.  Another  feature  of  Interest  is  the  fact 
that  as  longer  and  longer  life  times  are  required  of 
the  materials,  the  de -rating  in  maximum  allowable 
temperature,  for  a given  applied  stress,  is  much 
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SUMMARY 

Analytical  procedures,  developed  for  the  evalua- 
tion and  extrapolation  of  stress  rupture  data,  have 
been  applied  to  available  experimental  information  on 
several  high-temperature  structural  ceramics.  While 
it  is  evident  that  none  of  the  experimental  Information 
is  as  yet  sufficiently  extensive  to  fully  substantiate 
the  procedures  which  are  outlined,  the  results  are, 
however,  consistent  with  the  available  data.  Life- 
time estimates  resulting  from  these  procedures  indi- 
cate that  several  ceramic  materials,  recently  devel- 
oped, possess  high-temperature  stress  rupture  prop- 
erties that  may  be  of  significant  interest  in  the  design 
of  high -temperature,  long-lived,  load-bearing  ma- 
chinery. 
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Fig.  6(b)  10,000-Hour  stress  rupture  curves  for  sev- 
eral high- temperature  structural  ceramics  and  for 
an  early  (Niomic  80A)  and  an  advanced  (DS 
superalloy  as  determined  by  Eq.  (4)  and  available 
experimental  data. 

less  for  the  ceramics  than  for  the  alloys.  Also  of 
Interest  is  the  indication  that  sensible  strengths  of 
one  of  the  ceramics,  sintered  SiC,  are  predicted  for 
temperatures  well  in  excess  of  1600*C. 

For  applications  involving  inertial  forces,  de- 
signers generally  require  stress  rupture  data  re- 
ported in  terms  of  strength- to- weight  ratios  rather 
than  simply  strength.  Figure  7 corrects  all  the  ma- 
terials in  Fig.  6(b)  for  density  and  indicates  that 
much  of  the  disparity  in  low -temperature  performance 
between  alloys  and  structural  ceramics  is  removed. 

In  addition,  this  procedure  has  the  obvious  conse- 
quence of  moving  thecrossover  temperature  for  equiv- 
alent performance  to  much  lower  values. 
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